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Preface 


The aim of this book is to give a rapid introduction to differential 
equations which is elementary and which emphasizes the geometric 
nature of the subject. 

The book is elementary in the sense that it assumes only a know- 
ledge of the simplest properties of differentiation and integration. 
It is designed primarily for first-year college and university students 
who have not yet had courses in analysis or linear algebra. For this 
reason topics which are of only specialized interest (such as approxi- 
mation techniques and numerical methods of solution) and more 
advanced topics which require linear algebra and analysis (such as 
existence theorems, partial differential equations and series solutions) 
are either omitted or mentioned briefly in the guide to further 
reading. Occasionally, in order to give a precise statement of a 
definition or proof, it is necessary to use a property of continuous 
functions from an analysis course. The reader should accept such 
properties — which are usually ‘intuitively’ obvious — until he feels 
ready to look up further details in a book on analysis, such as A 
Preliminary Course in Analysis by R. M. F. Moss and G. T. Roberts 
(which has appeared in the same series, and is in the sequel referred 
to simply as Moss—Roberts). In this way the book should be suitable 
both for students specializing in mathematics and for students in 
other subjects who need differential equations as a tool. The latter 
group are warned, however, that they will not find the customary 
large number of examples of differential equations arising in biology, 
chemistry, economics, engineering and physics. The reason is that 
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such students are likely to study many examples in their own work, 
and are not likely to be helped by examples from a background 
which they do not understand. 

The approach is geometric, as can be seen from the large number 
of diagrams. The main reason is that most mistakes in differential 
equations result from blindly manipulating formulae without think- 
ing about the geometry involved. A geometric point of view is often 
necessary for applications, in which the required ‘solution’ of a 
differential equation is a picture rather than a formula. Finally, the 
geometric theory of differential equations has during the last hundred 
years developed from the elementary beginnings described in this 
book to an exciting area of research in which algebra, analysis and 
topology all play a vital role. 

I am grateful to Ronald Brown for his suggestion that such a 
book should appear in the Chapman and Hall Mathematics Series, 
to Oxford Illustrators Ltd for their help with the diagrams, and 
to the publishers and printers for their careful and sympathetic 
co-operation. 


R.L.E.S. 
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Symbol 


cos, sin 
cosh, sinh 


Glossary of 
Symbols 


Meaning 


cosine, sine 

hyperbolic cosine, hyperbolic sine 
exponential 

logarithm to base e = exp 1 

real part, imaginary part 

line of real numbers 

plane of pairs of real numbers 

set of infinitely differentiable functions f with 
fixed domain 

vector space 

domain of function 

region in R2 

belongs to 

function f which associates real number f(t) 
to each teD 

function f which associates value f(t) to each 
tin domain 


derivative of function 


partial derivatives of f(t, x) with respect to 
t,x 

positive square root (real number y >0 such 
that y? = x) of real number x > 0 


xi 
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“/x cube root (real number y such that y? = x) of 

real number x | 
el modulus (x if x >0 and — x if x <0) of real 

number x 30 


{x : P(x)} set of all x such that property P(x) holds 
(often context implies x € R) 8 


Introduction 


This book is based on a course of twenty lectures given to first-year 
mathematics students at the University of Warwick. One aim of 
the course was to provide an antidote to the ‘abstract’ courses on 
sets and groups which the students attended simultaneously in their 
first term. Another was to stimulate geometric intuition and to 
provide motivation for the courses on analysis and linear algebra 
which the students attended in the following term. 

These aims implied already that the course should be ‘concrete’ 
rather than ‘abstract’. There is, however, another reason for this 
approach. It is a common complaint among undergraduates that 
university mathematics seems unnecessarily abstract compared with 
the concrete mathematics they meet at school. Another aim of this 
book is therefore to deal with differential equations in a rather 
concrete manner, and to show how this leads naturally to such 
‘abstract’ concepts as continuous function, vector space and linear 
operator. 

In fact it is sometimes the school — rather than university — treat- 
ment of differential equations which is unnecessarily abstract. 
Consider for example the differentiation of functions like y = log x 
and y= sin x. It is familiar that the derivatives are 1/x and cos x. 
It may seem pedantic to point out that this is not quite accurate, 
since the function y=log x makes sense only for x>0, whereas 
its alleged derivative 1/x makes sense both for x > 0 and for x < 0. 

To see why such pedantry may be very necessary, consider the 
differentiation of a function of a function. The well-known chain 
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y= log x, x>0 


rule gives 
. (log sin x) = —— 
£ (log log x) = eae 
i HOB 1B SIO) Og 


However, these computations are somewhat abstract: an attempt 
to draw the graphs of the functions involved shows that the first 
two equations are often meaningless, while the third equation is 
always meaningless. Since y=log x makes sense only for x>0, 
the function y = log sin x makes sense only for points x lying strictly 
between an even multiple of m and the next odd multiple of 7z. 
On the other hand its alleged derivative makes sense for all points x 
(except for an infinite value at even multiples of 7). 


—_— == os a= oe o= ae as 


y = log sin x y = log log x 
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Similarly y = log log x makes sense only for x >1 whereas its 
alleged derivative makes sense for any x >0. Finally the function 
y = log log sin x has no graph at all because sin x < 1 and therefore 
log sin x <0. Whatever the formula involving y =log log sin x 
means, it is certainly not a piece of concrete mathematics. 

It may seem artificial to consider such a bizarre function, but 
the same kind of alarming situation can arise from the simplest 
differential equations. 

What is wrong with the following argument? ‘Since 
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the solution of the equation dy/dx = 1/x is y = log x.’ 
There are two mistakes. Firstly inspection shows that equally well 


d 1 
he (log (— ») nay 


so that there is another candidate y= log (— x) for the solution. 
Secondly either solution may be modified by a constant without 
affecting the derivative. Therefore any of the functions 


y=a+logx 
y=b+log(— x) 


(where a, b are real numbers) are solutions of the differential equation 
dy/dx = 1/x. This result clearly contradicts the statement, which is 
commonly made, that if (d/dx)( f(x)) = (d/dx)(g(x)) then there is a 
real number c such that f(x)=g(x)+ cc. The reason is that the 
statement is false unless more is said about the precise nature of 
the functions y = f(x) and y = g(x). It is not enough that they are 
represented by formulae; as we have seen in this area abstract 
formulae, however impressive, can be meaningless. 

Let us now get away from abstract manipulations of formulae, 
and consider the equation dy/dx = 1/x from a concrete point of view. 
The equation is an instruction to find functions y=/(x) whose 
graphs have slope 1/x at the point (x, f(x)). It is necessary to exclude 
the point x =0, since here 1/x is undefined and so the equation 
says nothing. At all other points the slope 1/x may be represented 
by a non-vertical direction. The diagram of slopes is symmetrical 
under reflection in the vertical axis. It is clear intuitively that the 
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region x > 0 can be filled by curves with the required slopes; these 
are in fact the graphs of the functions y = a + log x. A check shows 
that exactly one such graph passes through each point (x, y) with 
x>0. Similarly the curves filling the left-hand region x<0 are the 
graphs of the functions y=b+log (—x). Therefore the correct 
solution to the equation dy/dx = 1/x is 


y=a+logx, x>0 
y=b+log(— x), x<0 


where a, b are real numbers (sometimes called arbitrary constants; 
we prefer to avoid the use of this phrase because it is commonly 
asserted that if a differential equation involves dy/dx but no higher 
derivatives then its solution must contain only one arbitrary con- 
stant — the above example shows that this statement is false). The 
diagram shows that the differential equation dy/dx = 1/x does have 
solutions which fill up the two regions x >0 and x <0 in which 
the equation makes sense; however the formulae needed in the two 
regions are quite distinct. It is therefore necessary to specify a 
function, not only by a formula, but also by adding a precise state- 
ment of the values of x for which the formula holds. 
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y= b 
+ log (-x) 
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The aim of Chapter 1 is simply to introduce the minimum amount 
of language necessary to enable such precise statements to be made. 
In Chapter 2 a large number of examples are discussed from the 
commonsense viewpoint adopted above for the equation dy/dx = 1/x; 
in each example the method is first to draw slopes illustrating what 
the equation means, then to guess a formula for a solution in each 
region where the equation makes sense, and finally to check that 
there is a solution through each point of the region. The more 
formal discussion begins in Chapter 3. 

The reader who, as a matter of principle, reads only the intro- 
duction of any book, may like to amuse himself by considering the 
following examples, which I owe to Ronald Brown. Here and 
throughout the book, 1/x denotes the positive square root of a 
real number x > 0. 


1 What is wrong with the following examination question? 
Find the least value of the function y such that dy/dx = x — (1/x) 
and y=1 when x=1. 


2 It is often said that the general solution of the equation 
dy/dx = ./(1 — y?) is y=sin (a+x). This is false since for 
this it is not always true that dy/dx > 0. Make a correct state- 
ment. 
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2 
Why not attempt to evaluate | (sin ‘) dx by substituting 
—] j 


x= 1/u? 


Differentiate sin~! (2+ x?) and +/(1 — x?) + \/(x2— 1). Com- 
ment. 


An A-level question asks for the general solution of dy/dx = y/x. 
Show why it is not y=ax but 


y=ax, x>0 
y = bx, x <0 
where a, b are real numbers. 


Find a solution of dy/dx = 1 + 2(y/x) which has a minimum at 
x =1 and a maximum at x = — 3. 


CHAPTER ONE 


Differentiable 
functions 


We take for granted the existence of the real numbers, with all their 
usual properties. Individual real numbers may be denoted by symbols 
like 0, 1,2, 1933... Gy @ Oe By symbols like rf, x; ¥, x1, Xo, 222: 
We do not regard ‘infinity’ as a real number, nor do we allow 
symbols like t, x, y to denote complex values unless this is explicitly 
stated. 

In practice it is often preferable to think not in terms of algebraic 
formulae but in terms of geometric pictures. In order to do so, 
consider the set R of a// real numbers, and the set R2 of all pairs 
of real numbers. We think of R as an infinite line, and of a particular 
real number x1 as a point on the line. 


0 xy 
OTA Cs 


In the same way we think of R2 as a plane, and of a particular 


“4 
(0.y) if (x,y) y- coordinate UY 
hae (x,0) x- coordinate 
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pair (x1, yi) of real numbers as a point P on the plane. We call x1 
the x-coordinate and y, the y-coordinate of the point P. 

Some regions e” in R2 act as graphs determining yi uniquely in 
terms of x1. The condition for this is that e? should meet each 
vertical line in not more than one point. 


1.1 Definition 
A region e” in R2 is the graph of a function if no two points in of 
have the same x-coordinate. 


1.2 Examples 


We use two slightly different notations to specify regions in R2. 
For example ce” = {(x, y):x2 + y2 < 1} means that oe” consists of all 
points (x, y) whose coordinates satisfy the inequality x? + y? <1. 
Similarly 0” = {(x, y):x2+y2=1} means that e% consists of all 
points (x, y) whose coordinates satisfy x2 + y2 = 1. In this notation 
we write o% = {(x, y):some relation between x and y}. The other 
notation specifies both coordinates in terms of a parameter. Thus 
of = {(t?, t):t € R} means that oe” consists of all points (t?, t) where 
t is an arbitrary real number. In fact none of these three examples 
are graphs of functions, since in each case there are distinct points 
in ec? with the same x-coordinate. 


S={(x,y)ix?+y2 <1} LP={(x,y):x?+y2=1} S = {(t?, t):t € R} 


In each of the following examples, however, e” is the graph of a 
function. 

It is clear from these examples, and also from the discussion in 
the Introduction, that it is important to be precise about the values 
of x for which a given function makes sense. 


S = {(t,t):1 € R} SF ={(t,t):i—-l<t<tl} 


S = {(x, y)ix2+ y2=1 and y>0} SF = {(x, y)iy = x*} 


S = {(t3,t):t € R} S = {(x, y):y = 19% digit in 
decimal expansion of x} 
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The values of x for which the function is defined are said to 
form the domain of the function. This is the content of: 


1.3 Definition 


Let e¥ be the graph of a function. The domain Y of the function is 
the set of all real numbers which occur as x-coordinates. 


1.4 Examples 


oF = {(x, y):x2+ y2=1 and y > 0} is the graph of a function with 
domain —1<x< 1. Note that this graph could also be described in 
the form {(x, /(1 — x?)):-l<x< 1}. 

of = {(t3, t):t € R} is the graph of a function with domain R. It 
could also be described in the form {(x, x1/8):x € R}. 

These examples illustrate the fact that the graph of a function 
with domain J can always be described in the form 


oF = {(x, f(x)):x € BY. 


It is natural to think of f as a rule which associates to each point 
x € Da real number f(x), and to speak of the function f. 


“ 
‘< 


Z 


xt>f(x) Domain 9 
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We now have two ways of describing a function. As a graph: this 
leads to the common notation y=f(x), x €QZ. As a rule which 
associates a real number f(x) to each point x € Z: this leads to the 
less familiar notation xt+f(x), x € Z. The use of the special arrow 
++1s to avoid confusion with other contexts (e.g. — meaning ‘tends 
to’). These two descriptions are illustrated together: for some pur- 
poses the idea of a graph y=/(x) is more convenient, for other 
purposes the idea of a rule xt->/(x) is preferable. It will often be 
the case that f(x) is specified by an algebraic formula (for example 
I(x) = x*, f(x) = x), or that fis given a particular name (for example 
J (x) =sin x, f(x) = log x, f(x) = exp x). 


1.5 Examples 


A linear function has a graph of the form {(x, xb):x € R} for some 
real number b. We speak of the function y= xb or x+>xb. Note 
that under a linear function 0++0. An affine function has a graph 
of the form {(x, a+ xb):x € R} for some real numbers a, b. We 
speak of the function y=a+xb or x+>a+xb. If yo=a+xob the 
graph can also be written in the form {(x, yo + (x — x0)b):x € R} 
or y — yo = (x — xo)b. In both these examples the real number Bb is 
called the slope of the function. 


Y—Yo = (xX—X9)b 


(Xo,Yo) 


We recall that certain functions can be differentiated or integrated. 
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1.6 Differentiation 


Certain functions xt>/(x), x € D have a derivative xt>f'(x), x ED. 
The real number b = f’(xo) has the following significance. It is the 
slope of the affine function x++f(xo) + (x — x0)b which gives the 
closest approximation to x++>/f(x) near xo. If this approximation 
exists for all xo €QD then the function x+>f(x) is said to be 
differentiable. 


1.7 Remarks 


To make the above description of differentiation into a precise 
definition requires a precise explanation of the words ‘closest 
approximation’ and ‘near’ in 1.6. This requires the concept of a 
continuous function as given for instance in Moss—Roberts, Chapter 
2 (especially 2.1 and 2.4). The geometric idea behind the definition 
of differentiation is that if x+>/(x) is differentiable then each point 
(xo, f(xo)) of the graph has a non-vertical tangent. Every non-vertical 
line through (x0, f(xo)) is the graph of an affine function 


xt->f (x0) + (x — x0)b 


for some real number 6. Among these the tangent is the line giving 
the closest approximation to f(x). 


y = f(x) 


y =f(xXo) +(x—Xo)b 


Xo 


Note that, if the graph of a function xt->/(x) has a vertical tangent 
at (xo, f(xo)) for some xo € Y, then f is not differentiable. In this 
case the tangent is not the graph of an affine function. 
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differentiable not differentiable 


not differentiable not differentiable 


1.8 Integration 


A function x+> F(x), x € @ such that F’(x) = f(x) for all x € D is 
called an integral of the function f. If such an F exists we say that 
the function x+>f(x) or y =/(x) is integrable. 

To every method of denoting the function x:>/(x) there corre- 
sponds a notation for differentiation and integration. Thus we may 
write y= f(x) and dy/dx =f'(x). It is also common to write Df or 
df{dx for f’. A differential equation arises typically when a relationship 
is given between y and dy/dx, or equivalently, between f(x) and 
f(x). A function x++f(x) which satisfies this relationship is called 
a solution of the differential equation. 
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If the original equation is given only for certain values of x and y, 
for example for all points (x, y) in some region o%, then the graph 
of any solution will be contained in e”. Suppose that 7 is a region 
of the form {(x, y):x € D} and that the differential equation specifies 
f''(x)for all x € Z. Then any solution must have a domain contained 
in 2, and must moreover be differentiable, since otherwise f’(x) 
would not exist for all x €e Z. However, the fact that the function 
xt>f(x), x €Q is differentiable implies no statement whatsoever 
about the behaviour of functions given by similar formulae out- 
side Z. 

For example consider the equation dy/dx=y/x. Here f'(x) 
is specified only in the region {(x, y):x £0}, and the function 


y=ax, Be 3,8 
y = Dx, x<0 


is a solution. The equation x(dy/dx) = y is slightly different: now 
x =0 implies y = 0, so that there is a solution with domain R 


y= ax, x>0 
y=0, x=O0 
y = DX, hosel te 


This solution is not differentiable at x = 0 unless a = b; in accordance 
with the fact that the equation does not specify dy/dx along x =0. 
In both these examples it is not permissible to conclude from the 
equation that a= b. 


1.9 Example 


Consider the equivalent equations 


POse) of Sy, 
x 
The equation is given without restrictions on x, y; that is the relation 
holds for all (x, y)¢€R2. Any solution must be a differentiable 
function x -> f(x). A solution with graph containing the point 
(Xo, Yo) must have slope yy there. The required slopes can be illus- 
trated on a diagram. 

Clearly f(x) = 0 is one solution. It is plausible from the diagram, 
but not obvious, that there is a pattern of curves with the required 
slopes, and that each such curve is a solution. It is also plausible 
that there is a unique such solution through each point (Xo, yo). 
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1.10 Assumption 


We assume that the equation dy/dx =y has precisely one solution 
xt++f(x) through each point (xo, yo). 

This assumption implies the existence of certain differentiable 
functions, and can be proved using theorems on continuous functions, 
as for instance in Moss—Roberts, Chapter 4 (especially 4.5). For 
example the solution through (0, 1) has the graph {(x, exp x):x € R} 
and in fact this can be taken as the definition of the function 
xt> exp x. If on the other hand some other construction of exp 
is given, for example in terms of power series, then we could use 
the construction to prove assumption 1.10. The same graph through 
(0, 1) can also be written in the form {(log y, y):y > 0} and thus 
defines the function y+> log y. 

The usual properties of exp and log can be deduced from assump- 
tion 1.10 and the elementary properties of differentiation. Fix real 
numbers 5 and B, and compare the expressions 


y=exp (x +b) y= Bex x 
OF ante By ey ne 
dx dx 


Therefore y=exp (x+)) is a 
solution of dy/dx=y, and by 
1.10 the unique _ solution 
through (0, exp 5). 


Therefore ‘y= 8B exp x 1s. a 
solution of dy/dx = y, and by 
1.10 the unique — solution 
through (0, B). 


It follows that if B= exp b then B exp x = exp (6+ x) for all x ER. 
We deduce that exp (a+ b) = exp a.exp b for all a, b € R. Similarly 
log ab = log a+ log b for all a, b> 0. 

If n is a positive integer exp (n) =exp (1+ °°: +1) and hence 
exp (”) = e”, where e = exp (1). This suggests that we define e”, even 
when x is not an integer, by the equation e* = exp (x). 


1.11 Example 
Consider the equivalent equations 
P(x) oO) Jory dplax = 0, 


The equation is given without restrictions on x, y; that is the relation 
holds for all (x, y)é R2. Any solution must be a differentiable 
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function x!>/(x) with domain R. The required slopes can be 
illustrated on a diagram. 
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Clearly for any real number c there is a solution y =c. Similarly, 
if J denotes the interval a < x <b the equation 


FT teh 0, xel 


has solutions f(x) = c vith domain /. It is plausible from the diagram, 
but not obvious, that these are the only differentiable functions with 
domain / which satisfy the equation. 


1.12 Assumption 


We assume that if xt->f(x) is a differentiable function with domain 
a<x<bandf'(x) =0 for alla<x <b, then there is a real number c 
such that f is the constant function f(x) = c. 

This assumption can be proved using theorems on differentiable 
functions, as for instance in Moss—Roberts, Chapter 3 (especially 
3.4). It implies the uniqueness of certain differentiable functions. 
Suppose x+>/(x) and xt-> g(x) are two differentiable functions with 
domain a < x < bsuch that f’(x) = g’(x). Then the function x+-> h(x) 
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defined by /(x) = f(x) — g(%) has domain a < x < b and h'(x) = 0 for 
all a< x <b. Therefore by 1.12 there is a real number c such that 


I(x) = g(x) +c, a<x<b. 


Note that if the domain of x+>/f(x) and xt g(x) consisted of 
several disjoint intervals then the assumption would imply 


F(x) = a(x) + e(x) 


where x}->c(x) is a function which is constant on each interval. 
This is the situation in the example considered in the introduction, 
where the domain consists of two intervals x <0 and x>0, and 
any two solutions of the equation 


Hae) es, 
x 


differ by a function 
c(x) = @, x>; 
c(x) = db, x < 0, 


It is necessary to use such functions x+->c(x) in order to make 
correct the traditional rules about ‘arbitrary constants’ mentioned 
in the Introduction. 

In the sequel we make use without comment of the elementary 
rules for differentiation (in particular how to differentiate the sum 
of two functions, the product of two functions, and a function of a 
function) which should be familiar from any school treatment of the 
subject. These rules are set out explicitly for instance in Moss- 
Roberts, Chapter 2 (especially 2.5). The other important properties 
used in the sequel are those summarized in assumptions 1.10 and 
1.12. Any other property which is needed, will be mentioned as it. 
arises. 


Exercises 


1 Make use of the usual rules for differentiation to find derivatives 
for the functions +/x, log x, sin x, log 1/x, log log x, log sin x, 
log log sin x. 
In each case make a careful statement of the domain of the 
function, and the domain of its derivative. 
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Find integrals for the following functions with domains as 


indicated: 
iA 
(i) ie 


(i) * 


(ill) cos x 


sin x log sin. x 


x SO. 


Xie 0, 
all x, 


x #nz for integer n, 


Om ee Ton x oT. 


2nt < x < (2n + 1)z for integer n. 


Sketch the graphs of the functions 


Kite ey: xeER, 
xt (o/x)?, x20, 
xXt>V/(x?2), =x ER, 


Xl->—., x £0, 


(the four diagrams obtained should all be different). 


CHAPTER TWO 


Examples of 
differential 
equations 


Frequently the functions which occur in a differential equation are 
functions of time. For this reason we change notations slightly, and 
consider graphs {(t, f(t)):t € 2} of functions x =/(t) or tr~f(t) 
with domain Y. Suppose some relationship between f¢, f(t) and 
f'(t) is given; in what circumstances do we expect the existence of a 
function ¢+>/(t) satisfying the relationship, and how large is the 
domain of the function? It is instructive to consider more examples 
like 1.9. In the first few examples the equation specifies f’(t) in a 
certain region of the (rt, x)-plane. We adopt the foliowing procedure: 
firstly draw the lines of slope on a diagram, secondly sketch the 
curves with the given lines as tangents (the so-called solution curves), 
and thirdly use one of the curves to sketch the graph of a typical 
solution of the equation. 

In many cases the shape of the curves suggests immediately how 
to guess the equation x = f(t) of a solution. Once a guess has been 
made it is usually easy to check whether the function x = /(f) is in 
fact a solution, and whether there is a solution through each point 
(t, x). The discussion of explicit methods for deducing the equation 
x = f(t) of a solution is postponed until Chapter 3. 


2.1 Example 
f’(t)=1, also written dx/dt =1. 


The equation asserts that the slope at each point is 1. The diagram 
suggests that every solution has the form x=f+c, where c is a 
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real number, and that there is a unique solution with domain R 
through each point. In fact both these statements follow from 
assumption 1.12. 


Slope lines, solution curves and typical solution of f(t)=1 


2.2 Example 
S(t) =2t, also written dx/dt = 2t. 


In this case there is again a unique solution with domain R through 
each point. Clearly x={f? satisfies the equation; therefore by 
assumption 1.12 every solution with domain R has the form 
x =1t2-+c where c is a real number, and there is one such curve 
through each point (f, x). 


2.3 Examples 

Consider the differential equation 
I @=1, t>0 
f@=-1, t<0. 


In this case f’(1) is undefined at t=0, so that we expect solutions 
with domain R* = {t:t 40}. In the region t>0 there is a unique 
solution x=t+a, t>0O through each point; in the region t<0 
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there is a unique solution x = —t+56, t<0 through each point. 
However, a typical solution with domain R* has the form 
x=t+a, t>0 


x= -—-t+5, t<0 


where a, 6 are real numbers. Again these statements follow from 
assumption 1.12. 


The situation is similar to the equation f’(t) = 1/t, t 4 0 discussed 

in the introduction. Here a typical solution has the form 
x =log t+a, t>0O 
x=log (-1) +), t<0 

where a, b are real numbers. 

Notice that the equation f’(t) = t/t, t 4 0 is not the same as example 
2.1, since the slope is undefined at t=0. A typical solution has the 
form 

x=t+a, t>0 

x=t+5, t<0. 
This distinction can be important in applications, where an expres- 
sion like t/t occurs in the middle of a calculation and is cancelled 
unwittingly. 


2.4 Example 
ft) =f) —2, or dx/at =x — 2: 
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The situation is almost the same as that in 1.9. By assumption 1.10 
every solution has the form x =2+ ce! where c is a real number, 
and there is a unique solution through each point (f, x). 


2.5 Example 
fOH=2V f(t), or dx/dt =2/x. 


Here, and in the sequel, »/ denotes positive square root. The equation 
implies that x >0. In the region x >0 there is a unique solution 
through each point. The function x =(t—c)?, t>c satisfies the 
equation, and has domain t > c for some real number c. However, 
in the region x > 0 there are solutions of the form 

x =(t—c)?, t>c 

x=0, t<sc 
with domain R, as well as the solution x =0. Note that through 
any point on the line x = 0 there pass an infinite number of solutions. 
The reader may suspect, correctly, that this is connected with the 
fact that at x =O the function xt 2,/x is not differentiable. 


2.6 Example 
t+f(t)f'(t) =0, or t+ x(dx/dt) =90 


x=- ¥(b-12) 


In this case the slope is undefined when x = 0. The diagram suggests 
that curves t? + x2 =c will satisfy the equation for any real number 
c > 0. This is in fact the case because (d/dt)(t? + x?) = 2[t + x(dx/dt)]. 
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Thus typical solutions are x = »/(a— t?) with domain {t:t? <a}, 
x=0 with domain {t:t=0}, and x= —+/(6—1#?) with domain 
{t:t2 <b}, where a, b are positive real numbers. There is a unique 
solution through each point (t, x) with x #0. 


2.7 Example 
tf'(t) + f(t) =0, or x +t(dx/dt) = 0. 


nt A NA, oe Wi * ne 
ie, tA 
Mi, \ f ae ag 
\ b>0 a<0 


The slope is undefined when ¢ = 0. The diagram suggests that curves 
tx = c will satisfy the equation for any real number c. This is in fact 
the case because (d/dt)tx = x + t(dx/dt). Through each point of the 
region ¢ < 0 there is a unique solution of the form x = b/t with domain 
t<0. Typical solutions with domain R have the form 


a 

x ae. t>0 
t 

x= 0, t=0 
b 


where a, 5 are real numbers. 


2.8 Example 
f(t) =3f(), or t(dx/dt) = 3x. 
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The slope is undefined when ¢ = 0. The diagram suggests that curves 
x = ct? will satisfy the equation for any real number c. This is in 
fact the case because t(d/dt)(ct®) = 3ct?. Through each point of the 
region t<0 there is a unique solution of the form x= bt? with 
domain ¢t <0. Typical solutions with domain R have the form 

ye ar, t>0 

eee 0, t= 

or. t<0 
where a, 6 are real numbers. 


2.9 Example 

fO=(2t+/Oy/[t-—2f(], or dx/dt = (2t + x)/(t — 2x). 

The slope is undefined along t=2x, zero along x= —2t, and 
1 along t= — 3x. Other slopes are 2 along x =0, and —4 along 
t =0. The equations (d/dt) (x? + t?)=2x(dx/dt) + 2t=2(x? + 2?)/(t-—2x) 
imply that as ¢t increases the distance from (t,x) to (0,0) increases 


x=—2t nae 


Slope lines and solution curves for at gas 
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or decreases according as t>2x or t<2x. Therefore the curves spiral 
in to the origin in a clockwise direction as shown. 

Typical solutions are given by fragments of spirals chosen so as 
to be graphs of functions (as in definition 1.1). In this case we make 
no attempt to guess a formula for x in terms of t. The reader may 
suspect, correctly, that to do so conveniently requires polar co- 
ordinates or complex numbers. 


2.10 Example 
t (2x3 — £3)(dx/dt) + x(2t? — x°) =0. 
P20 


sign of slope slope lines 


CIS 


solution curves 
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The slope is undefined along t=0 and t= x4/2, zero along x =0 
and x= 14/2, and —1 along x =¢t and x= —t. 


2.11 Examples 


Examples 2.9 and 2.10 illustrate the fact that a sketch of the solution 
curves may well give useful and detailed information, even though 
it is not possible easily to write down formulae for the solution. 
It might be thought that every differential equation decomposes 
regions of the (t, x)-plane into curves which hardly ever cross each 
other. This is true for sufficiently simple equations, but not for the 
more complicated kinds of equation which occur in the following 
examples. 

Consider the equation (dx/dt)? = t. It implies first that t > 0, and 
secondly that at the point (ft, x) there are two possible slopes 4/t 
and —4/t. The corresponding diagram of solution curves shows 
two curves through each point (t, x), t > 0. Notice that this pattern 


x 


can be got by superimposing the solution curves for the equation 
dx/dt = 4/t and the equation dx/dt = — +/t. 

A similar situation arises from the equation (dx/dt)? = 4x. The 
equation implies that x>0 and dx/dt= +2+/x. In the region 
x>0O the solutions to dx/dt=2,/x (compare example 2.5) and 
dx/dt = —2/x are superimposed. 

Differential equations involving (dx/dt)? but no higher powers 
are said to be of the second degree. Such equations frequently 
determine a region of the (t, x)-plane in which there are no solutions, 
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x 


{ 
Wind 


t 


and a region in which there are precisely two solutions through 
each point (obtained by superimposing solutions of two equations 
of first degree). 


2.12 Examples 


Consider the equation (d2x/dt?) =0, or f’(t) =0. 


This equation involves the second derivative of t'->/(t), and is 
therefore said to be of second order. It is satisfied by any affine 
function 


x=a+tb 


where a, b are real numbers. There is an infinite collection of solution 
curves through each point (ft, x), since every non-vertical line in the 
(t, x)-plane is the graph of a solution, and hence there is no con- 
venient diagram of solution curves. 

A better approach to this equation is to substitute y = dx/dt, and 
to write the equation 


Now a solution can be regarded as a pair of functions t+ x(t), y(t) 
and represented by a curve in the (x, y)-plane. Clearly the solution 
curves have the form y=a; the equation dx/dt = y shows that x 
is increasing if a>0O, constant if a=0, and decreasing if a<0. 
Solutions are illustrated in the diagram overleaf: points on opposite 
sides of the line of rest points move in opposite directions. 
This method of replacing a second-order equation by two first- 
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J ; 


— ee emnnieenemomnmmmmmnntind 
—etitidneeiniemmned 


order equations involving t, x, y, dx/dt, dy/dt is adopted in Chapter 
4. In Chapter 3 we restrict attention to first-order equations. 


Exercises 
1 Sketch solution curves of the differential equations 


dx 


eee Dt 
aR 
dt 

dx 
pare 
aie 


2 Let |+| denote the modulus of ¢, and 4/ positive square root. 
Explain the differences between the equations 


dx \ 2 dx\2 dx dx dx 
(7) =, (<) =|t], aay ag mats ap wa 
and sketch solution curves for each (the five diagrams obtained 
should all be different). 


3 Consider the differential equation 


Prove that there are no solutions through points in the region 
{(t, x); x2 < 41}. What can be said about the solution curves 
in the rest of the (t, x)-plane? 


CHAPTER THREE 


Solution of 
first-order 
equations 


The method used to solve the differential equations of the previous 
chapter was essentially to guess the formula for a solution and then 
to check that it works. Often however it is possible to deduce the 
form of a solution by using elementary integration. This is particu- 
larly the case for equations of the first order, that is, involving 
t, f(t) but no higher derivatives. We make use of the usual rules for 
integration (including substitution and integration by parts) as well 
as the standard properties of trigonometric and exponential func- 
tions. A careful discussion of all these matters is given in Moss— 
Roberts, Chapter 4 (especially 4.4, 4.5, 4.6). 


3.1 Equations of the form dx/dt = g(t) 


Examples 2.1-2.3 are of this form, in which dx/dt depends only 
on ¢ and not on x. The examples illustrate that the domain of the 
function g may not be the whole line R, and that to obtain a solution 
given by a single formula it is necessary to restrict to an interval 
p<t<gq in which g is well behaved. 


We adopt the convention, when referring to an interval p < t <q, 
that the interval may in fact be unbounded at one end (p<? or 
t<q) or at both ends (the whole line R); however, we prefer not 
to use picturesque symbols like ‘gq = + 00’ to describe this situation. 
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Theorem. Let g be an integrable function with domain p<t<q. 
Then the differential equation 
dx 
eee | ee JOE 
Fi 
has a unique solution x = f(t) through each point (to, Xo) with p < to <q. 
The domain of f is the interval p< t <q. 


Proof. Since g is integrable there is, by 1.8, a function G with 
domain p<t<q such that G’(t)= g(t). Let c= xo — G(to). Then 
(t) = c+ G(t) gives the required solution. If A(t) is another solution 
with A(to)= xo, and k(t)=/(t)—A(t), then k is a function with 
domain p <t<q which satisfies the equations k’(t)=0, k(to) =0. 
By assumption 1.12 this implies k(t)=0. Therefore the solution 
x = f(t) is unique. 

Remark. The lack of content in this proof is due to the fact that 
we defined integration in 1.8 as the opposite of differentiation. If 
integration were defined in some other way - for example by cal- 
culating areas or by integrating power series term by term — the 
proof would be more difficult and the theorem more significant. 
The application of this theorem depends on being able to compute 
the integral G in a specific case; we give some examples of this, 
using elementary facts about integration. 


3.2 Examples 
Let a, b be real numbers, and consider the equation 


dx 
=at+b. 
Wy a 
By assumption 1.12 there is at most one function ¢}--> G(t) with 
G'(t)=at+b and G(0)=0, so clearly it must be t+>4at? + bt. 
Therefore every solution to the equation has the form 


x = }at?+ bt+c 


where c is a real number. 
Now consider the equation 


dx 1 


Ua ~l<t< I, 
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By assumption 1.12 there is at most one function t+>G(t) with 
G'(t) = 1/+/(1 — t?) and G(0) = 0. We use the standard notation 


G(t) = | Par; Waa du 


for this function. Since —1<t<1 we can write u=sin @ with 
— }n <0 < 4m, and obtain 


sin-'¢ 


t 
l | Ae 
| se u-] ery dil donl an 1 ¢, 
0 0 


Here sin-1 ¢ denotes the unique number @ such that sin @= tf and 
— 4 <6< 4m. Thus every solution to the equation has the form 
x=c+sin-!f, —l<t<1 
where c is a real number. 
Consider next the rather similar equation 

nl RN 

dt /(t?—1) 
Here the domain of the function g(t) = 1/+/(t? — 1) consists of two 
distinct intervals which must be considered separately. Define 


t<-—lort>l. 


cosh v = $(e” + e-*) 

sinh v = $(e” — e~”). 
The properties of exp in 1.9 and 1.10 imply that (d/dv)(cosh v) = sinh v 
and (d/dv\(sinh v)=cosh v. Moreover (cosh v)?—1=(sinh »v)?. 
If u>1 there is a unique real number v >0 such that cosh v=u, 
and we write v = cosh! wu. These facts suggest the following substi- 
tutions. 

(i) In the interval u> 1 let «=cosh v, v>0; then 


cosh~' t 


t 
1 a mat wl 
| erp tia | d= comm t 
1 0 


and so every solution has the form 


x=a+cosh-!t, i> i; 
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(ii) In the interval u< —1 let —u=cosh v, v>0; then 


cosh-! (—t) 


t 
I : 
J LRT) al Lice 1(—f) 


and so every solution has the form 


x = b—cosh-! (— 4), t<—l1. 


| 


si 


SS 
| 
<= se Ss 


Goran 


Solutions to Solutions to 
pI ue 
dt /(l—?t?) dt /(t?—1) 


Consider the equation dx/dt =1/t, t#0. Again the domain of 
g(t) = 1/t consists of two distinct intervals t>0 and ¢<0. By 1.10 
every solution in the region ¢ > 0 has the form x = a + log ¢. Similarly 
every solution in the region t <0 has the form x = b + log (— 1). 

The substitution ¢t — qg for t, shows immediately that every solution 
of the equation dx/dt = 1/(t — q) has the form 


x=a+log (t—q), gat, 
x=b+ log (q—12), t<q. 
Finally consider the equation 


dt (p—tq-ty 
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The domain of g(t)=1/(p—t)(q—t) consists of three intervals 
t<p, p<t<q, and q<t. Since 


oe! (*- -*5] 
CP ERG Eee Vee | oe 
every solution has the form 


1 
x= ¢+-—— [log (q—¢) — log (p—# 
rae ip t)] 
my one log (=) for f <P, 
q—p | doa 4 
1 
has [log (q — t) — log (t — p)] 


] q—t 
Btu log (=) forp<t<4q, 
q—p a0 


b at [log (¢ —q) — log (t — p)] 


Byes log (‘= 4) for g <t. 
q—Pp ip 


3.3 Equations of the form dx/dt = h(x) 


Examples 2.4-2.5 are of this form, in which dx/dt depends only 
on x and not on ¢. The examples illustrate that the domain of h 
may not be the whole line R, and that to obtain solutions given 
by a single formula it is necessary to restrict to an intervala<x<b 
in which h is well behaved. Moreover it may be necessary to restrict 
to an even smaller interval to make sure that there is a unique 
solution through each point (fo, Xo). 

We first observe that if h(c)=0 then x=c is a solution. This 
makes it reasonable to restrict attention to intervals a<x<b on 
which A(x) > 0 or h(x) < 0. 


Theorem. Let h be a function of constant sign with domaina<x <b 
(that is, either h(x) > 0 for alla< x <b, or h(x) <0 for alla<x <b). 
If 1/h is integrable, the differential equation 


has a unique solution x = f(t) through each point (to, Xo) witha < xo < b. 
The domain of f is some interval containing to. 
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Proof. Note first that by 3.1 the equation 


wi oa a<x<b 

dx: Wx) 

has a unique solution ¢t = k(x) through (to, xo). Thus k is a function 
with domain a< x <b, and k'(x) is of constant sign. We now use 
the inverse function theorem proved in Moss—Roberts, Chapter 3 
(especially 3.5) to assert the existence of a unique function f, with 
domain an interval containing fo, such that, for alla<x <b, 


(i) x = f(t) if and only if t = k(x), and 


(1) x =f(t—c) satisfies dx/dt = h(x) if and only if t=c+ k(x) 
satisfies dt/dx = 1/h(x). 


We conclude that every solution in the region a<x<b has the 
form 


x=f(t—c). 


Note that in practice the formula for the inverse function fis usually 
immediate, so that the full strength of the inverse function theorem 
is not needed. 


3.4 Examples 


Consider the equation 


A comparison with the proof of 3.3 shows that 
oa) =X 
h(x) 

and therefore solutions of dt/dx =x in the regions x >0 or x <0 


have the form 
t=c+4x%, 


Therefore through each point (fo, xo) with xo> 0 there is a solution 
of the form 
x = o(2t — 2c) 


with domain {t € R:t> c}. 
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Consider next the equation dx/dt = }(x? — 1). There are solutions 


x =1and x= —1 because for these values h(x) = $(x2 — 1) 1s zero. 
A comparison with the proof of 3.3 shows that 
1 P l 1 


Wx) (x— 1)\(x + i i eee 
and that therefore the solutions of dt/dx = 1/h(x) in the regions 
x<—1orx>1 have the form 


t=c+ log esti 
x+1 
This equation is equivalent to e’*=(x—1)/(x+1), so that 
1 + ef-¢ = x(1 — et-¢). Through each point (to, Xo) with xo > 1 there is 
a solution of the form 


I Nik MI Ws ak 
ee. ee ete 1 — aet 
where c is a real number and e~¢=a> 0. Note that the domain of 
this function is {t €e R:t< — log a}. 
A similar calculation in the region —1<x<1 shows that 
t=c+log (1 —x)/(1 +), and hence 


MG? oust aiVADN Sigel’ ad 
ee ee ee 
is a solution with domain R for any a> 0. 
The solutions to this equation are illustrated in 4.1. 
Now consider the equation dx/dt=./(l— x?) —l<x<l. 
There are solutions x =1 and x= —1 because for these values 


h(x) = »/(1— x?) is zero. A comparison with the proof of 3.3 
shows that 


Oe 
h(x) VU — x?) 


and that therefore the solutions of dt/dx =1/h(x) in the region 
—1<x<1 have the form 


t=c+sin—! x, —l<x<l. 
This equation is equivalent to 


x =sin (t—c), c—tr<t<c+ 3, 
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and yields solutions with domain c—47<t<c+4m. Through 
each point (fo, xo) with —1<xo<1 there is a unique solution. 
However, through points (to, — 1) and (to, 1) there are many solu- 
tions, for example every real number c determines a solution 


x=-—1 t<c— $n, 
x = sin (t —c) c—dr<t<ct na, 
x=1 ctim<t. 


For a further example consider the equation dx/dt = 3x2/3, Here 
x?/3 denotes the unique cube root of the real number x2. There is 
a solution x =0 because for this value x2/3 is zero. Each point 
(to, Xo) in the region x >0 lies on a unique solution of the form 
x = (t—q)* with domain t>q. Each point in the region x <0 lies 
on a unique solution of the form x =(t—p)? with domain ¢ <p. 
However if p <q there are solutions with domain R of the form 


a 4 ne 


egy 
and so uniqueness fails for points on x = 0. 


3.5 Equations reducible to the form 3.1 or 3.3 


The first-order equations considered in: 3.1 and 3.3 are extremely 
special. A more typical equation would have the form 


dx 
dt =f(t, x), (t, x) ee 
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where e” is a region in the plane R? and (ft, x)t~>/(¢, x) denotes a 
function which associates, to each point (t, x) of e”, a point f(t, x) 
of R. Occasionally, however, it is possible by cunning manipulation 
to reduce the above equation to one of the form 3.1 (in which f(f, x) 
does not depend on x) or of the form 3.3 (in which f(t, x) does not 
depend on t). We give some examples. 

Consider an equation of the form 


dx 

ai oe OED 

= g(O)h(x) 

where g is a function with domain p <t<gq, and h is a function 


with domain a < x < b. If there are solutions for the equations 


ep p<t<q, as in 3.1 
dt 

dx 

a AX), A OLAS In 3.5 
dy 


then x =x(y(t)) is a solution of the original equation inside the 
rectangle p<t<q,a<x<b. For instance the equation 


dx 
a 2tx 


can be solved in this way. The equations 


have solutions y=?t?2+c and x=be¥ respectively. Therefore a 
solution to the original equation dx/dt = 2tx is 


x = bet*+e = (be*)et* = aet" 


where aE R. Equations of this type are sometimes said to have 
variables separable. 


Consider now an equation of the form 


dx 
a = g(u) 
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where u=x/t and g is a function with domain a<u<b. This 
equation is therefore defined for all points (t, x) in the region strictly 
between the lines x = at and x = bt. Note that (0, 0) is excluded, so 
that the region consists of two distinct parts. If x = ut we have 


dx du 


g(u) wy ag t+u 
and hence an equation 
du 
FF aay (g(u) — u) 


Zi 


Solutions to 
dx _x?+?? 
dt t? 


which is of the ‘variables separable’ type (note that t 40, so that 
1/t is well defined, throughout the region in question). For example 
the equation 


is of this type, since x = ut gives 


dye Dep an yt | 
dt 


— =— (u—1)?. 


The equation dy/dt = 1/2t has solutions of the form 
y= log|t|+ce 
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in the regions t>0 and t<0. The equation du/dy =(u— 1)? has 
a solution u = 1, as well as solutions u=1—1/(y+q). Every point 
(xo, to.) with to #0 lies on a solution of the original equation 
given by a formula of the following type (we set a= 2(c+ q)). 


2t 


If xo > to >.0, x = t——___-—— with domain 0 < t < e~4. 
a+logt 

Ifxo=to>0,x=t with domain 0 < t. 

Iffo> Xo>Ox=f TON ON with domain e~@ < t¢. 
a+logt 

2t ' ; 

If xo < to < 0, x = t -——_—_—._ with domain —e-* <¢ <0. 
a+ log (— 1?) 

If x9 = to <0, x= with domain t < 0. 

If to << x0 <0, x= EN oN with domain t < — e~@. 
a+log(—t) 


The reader should use these formulae to make an accurate version 
of the rough sketch of solutions on the opposite page. Equations 
of this type are said to be homogeneous, but this is only one of the 
many distinct meanings of the word ‘homogeneous’ in connection 
with differential equations. 


3.6 Exact equations 


As in 3.5 we consider a function (t, x)+~/f(t, x) which associates, 
to each point (t, x) €0%, a point f(t, x)¢R. We use without 
proof the fact that certain such functions are differentiable, and 


yield ‘partial derivatives’ bd and e Thus tt» = us (t, Xo) 


is defined to be the derivative of t+->/f(t, xo), for fixed xo. yi 
Xb> 2 (to, x) is the derivative of xt~>/f(to, x), for fixed fo. 


If x =g(t) then f(t, x) can be regarded as a function of ¢ with 
derivative 


d 2+ of dx 
it CME AD) at a 
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This fact can often be exploited as follows. Consider a differential 
equation of the form 


We NB Be 
dt 
If it is possible to find a function (t, x)t»f(t, x) such that 
(t, x)=a(t, x) and .: (t, x) = b(t, x) then the equation can be 


written 


d 
5 L62)) = MO) 


and solved by the method of 3.1. We give two examples. Consider 
the equation 


d. 
(pt + qx) + (qt + sx) “, =0 
where p, g, s are real numbers. It can be written in the form 
© pt + gtx +} sx) =0 


and therefore there are solution curves of the form 
pt? + 2qtx + sx*=c 


for any real number c. Fragments of these curves are graphs of 
solutions x = g(t) of the original equation. 
Note that an attempt to apply the same method to the equation 


(pt + qx) + (rt + sx) <0 
will fail if r#~q. This illustrates the fact that it is only in certain 
lucky situations that the above method can be used. 
For another example consider the equation 


dx : 
2x it + x2 cos t= 2t esin Ff, 


It can be written in the form 


dx 
x2 e8in t cos ¢ + 2xesin t a == 2t 


. (x2 esin t) = 27 
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and therefore there are solution curves of the form 

Mee GP cpa on 
for any real number c. The fragments of these curves given by 


x=et sin ¢ ./(t2 + a), t2> — a, 
x= —e teint ./(t2+5), 2-5, 


are solutions to the original equation. The method used in this 
example consisted of multiplying through by a non-zero ‘integrating 
factor’ e®i2 4, designed to make the equation exact. There is one 
class of differential equations for which this method always works; 
these are the so-called Jinear equations. 


3.7 Linear equations 
A first-order equation is Jinear if it can be written in the form 


* + x(t) = ht) 

where g, A are functions with domain p < t <q. We wish to investi- 
gate under what circumstances there exists a unique solution through 
each point (to, Xo) with p < to <q. 

We first assume g integrable. Then there is a unique function G 
with domain p<t<q such that G’(t)=g(t) and G(to) =0. Since 
e@) is non-zero for all p<t<q, the differential equation is not 
changed by multiplying through by the integrating factor e¢. 
Therefore the original equation is equivalent to 


x eG g(t) + eG 3 = eG h(t) 


© (e010) = h(t) 80. 


We now assume he integrable. Then there is a unique function 
k with domain p <t<q such that k’(t)=h(t)e©® and k(to) =0. 
By 3.1 there is a unique solution through (fo, xo) of the form 

x eG = x9 + k(t) 
x = e~& (xo + k(t)). 


The solution is a function with domain p <t <q. 
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Remarks. The usual notation for G and k in terms of integrals is 


t 
G)=] eu) du, 


to 


k(t)=| e&@ h(u) du. 


ot em, nw 


The most important special case is when g is constant. Then the 
equation has the form 


dx 
on + xa = h(t) 


for some real number a. Every solution can be written in the form 


t 
x = e74 (¢ + i e@ h(u) du 
to 


where c is a real number; this solution passes through (to, xo) if 
c = xpe™o, 

The theory of 3.1, 3.3, 3.5 and 3.7 may raise the following 
question. Clearly it is of interest to know whether a given differential 
equation has a solution through a given point (to, xo). But why 
is it of interest to know whether or not the solution is unique? 
There are two main reasons. Firstly, in applications where x is 
a variable depending on time f¢, it is important to know whether 
or not the differential equation and the value xo at time fo deter- 
mine x uniquely for all t>to. Secondly, differential equations are 
often solved in practice by programming a computer to make 
successive approximations to the solution; if the solution is not 
unique the computer may have a nervous breakdown trying to 
find it. 


3.8 Equations reducible to linear form 


Occasionally an equation which is not linear can be reduced to a 
linear equation by a suitable trick. For example consider the 
‘Bernoulli equation’ 

dx 


sit = xn 
i +x g(t) = x" h(t), x #0. 
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Let y= x1-", which makes sense because x #0. Then 
in & 
a. (1 —n)x-” oa (1 — n)A(t) — (1 — n)x!-* g(t) 
so that y is a solution of the differential equation 
d 
=~ a — Dy gt) = — (2 — DAC) 


which is linear. 
For another example consider the equation 


t 4 +xe(th=h(t), t>0. 


d 
Let t = e“, which makes sense because ¢ > 0. Then 
dx _dxdt dx , 
du dat du dt 


so that x is a solution of the differential equation 


a +. x g(et) = h(et) 
du 
which 1s linear. 

In both these examples the substitution which makes possible 
the reduction to a linear equation depends crucially on the restric- 
tions imposed (x ¥ 0 in the first example; ¢ > 0 in the second). 


3.9 Examples 
Consider the equation 


is an xo 0 
Axe 
This arises in economics, where c is the total cost of an output x, 
dc/dx is the so-called marginal cost, and c/x is the so-called average 
cost. This equation is homogeneous. The substitution c= ux gives 
du/dx =0 and hence every solution has the form u=constant. It 
follows that if the marginal cost equals average cost at every output 
then the average cost is constant. 
Consider the equation L(dx/dt) + xR = v(t). This arises from an 
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electric circuit with inductance L > 0, resistance R > 0, and voltage 
v(t) at time t. The solution x(t) is then the current at time ¢. This 
equation is linear. The method of 3.7 gives, writing a = R/L, 


dx 1 
es + xa =F v(t), 


a. 1 
rat at) — — (ft) eat 
Fe eed az veem, 
t 
xet%#=c+ : | v(u) e@ du. 
0 


Suppose for example that v(t) is a constant input E started at time 
t = 0 and switched off at time ¢ = t1. There are then two cases. 


For 0 <¢t< t; choose the solution with x =0 at t=0. This cor- 
responds to choosing c =0 and v(u) = E. Therefore 


t 
E E 
at — au di me TOM IU OG i: 1 
xe I { ert du 7 (e ) 
0 


and hence 
E : 
4 ee Rt/L) for O<t< fy. 
For t1 < t choose the solution with x = (E/R)(1 — e~%1) at t= 4%. 
This corresponds to choosing c=(E/R)(e%:—1) and v(u)=0. 


Therefore 


E 
x= pa a ee lehte. for nn <4, 
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Exercises 


1 Find every solution t++ x(t) of the following differential equa- 
tions, as well as the solution through (1, 0). 


i) i a eal 
Hh + 4f, (il) = bet, beR, 
1 inl ee l 
id) TR OY) ae ey 
dx a ee NE 


2 For each of the following differential equations find every 
solution t+>x(?t), and state its domain. Is the solution through 
(0, 1) unique? 


e dx 2 ee dx x 
array Ph sla R 
(i) 7 x2 — 3x + 2, (11) it bet, bER, 
A OR A Et ai 
(iii) a ieee)", (iv) ar Mie, ee hy 


(v) om 2s, eG, (vi) — oe ta x #(n+ 4)z. 
3 Find all solutions of the differential equations 
dx 
(i) (3t2x — tx) + (213x?2 + t3x4) a 0, 
a dx 
(ii) (1 + 2x) + (4-1?) ita ie a 


x 
(ini) - & = cos 7m i >'0, 


(iv) (t? — x?) = a ft i> x> 0, 


(v) et & + 3xe3t = 2t, 


beh (ee Sey Ete . 2 
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Solve the differential equations 


i ae nae! 
(i) Fase 


(ii) © + xtant=0, a attes 

URED le ved it) ets Mee 
dt f 2 vy 

dx 

(iv) ¢3 Gee ee ae t>0, 
dx F 

(v) G+ 2tx + tx = (), x iD> 0; 

(vi) 1 + xlog t= 22, t>O0. 


Show how to solve an equation of the form 
dx 


de 


+ x g(t) = h(t), t<0O 
by a substitution t= — e%, 


Show that x =7¢? and x=74 are solutions of the differential 
equation 


d2x dx 

te __ — 6¢ —- + 12x = 0, 
Maa? Hida 

Deduce that if p, g, r, s are real numbers, and 


ptie+qtt t>0 
f= 


rt3 + st4 t<0 


then x = f(t) is a solution. What restrictions are imposed on 
DP, q, r, 8 if it is required that f be 

(i) twice differentiable, 

(ii) three times differentiable, 

(iii) four times differentiable? 


CHAPTER FOUR 


: Autonomous 
systems 


In the previous two chapters solutions to a differential equation 
dx 
ed ae , — 
aed) 


were considered as functions f+ x(t), or as curves in the (f, x)-plane. 
This description has certain disadvantages: 


(i) unless the function f has an especially simple form it may not 
be possible to find explicit formulae for a solution; 


(ii) such exact information may in any case not be required for 
the problem under discussion, and it might suffice to give a simpler 
description; 


(iii) a physical problem might involve two functions x, y related 
by differential equations 


d 
7 = g(t, x, y) 


nae 
dt raat A(t, x; y) 


so that solutions must now be represented by curves in the three 
dimensional space with coordinates (ft, x, y). 

The equations in (iii) may be imagined as controlling the motion 
of a particle in the (x, y)-plane. The velocity of the particle, if it 
reaches (x, y) at time f, has components dx/dt, dy/dt. In some physical 
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problems the velocity depends only on (x, y), not on the time ¢; 
that is, the direction and speed of travel depend only on the position 
(x, y) of the particle, not on the time at which it reaches that position. 
This situation corresponds to a pair of differential equations 


dx 
ite 2(x, y) 
dy 
dt HAG h(x, y) 


where g, / are functions defined in a region e” of R2. Such a pair 
is called an autonomous system in two variables. A solution of the 
autonomous system is a pair of functions t+->(x(t), y(t)). In the 
same way, an equation 
dx 

Fe 
of the type considered in 3.3, is called an autonomous equation 
(or an autonomous system in one variable). The essential feature, 
in both cases, is that the right-hand side does not depend on ¢. 

The frequent occurrence of autonomous systems in two variables 
leads to a search for a simpler method of representing geometrically 
their solutions. First, we consider an example of the corresponding 
situation for one variable. 


4.1 Examples 

Consider the equation 
dx F 
a = oat tS 
= i(x?—1) 


It is autonomous, because the right-hand side is independent of 1, 
and explicit solutions were obtained in 3.4. For many purposes, 
however, the following analysis would be sufficient. 

At x= +1 we have dx/dt=0, so x=1 and x= —1 are rest 
points; for x >1 or x < —1 we have dx/dt>0, so at these points 
x increases with time; for —1<x<1 we have dx/dt <0, so at 
these points x decreases with time. The left-hand side of the diagram 
represents the information given by the explicit formula: x =1, 

= —lor 
ol ae 


Sygate 


x 
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obtained in 3.4. The right-hand side of the diagram represents 
the information obtained in the cruder analysis above. Both diagrams 
display equally the most important facts about the equation, namely 
that x=1 is an unstable solution (points near x =1 move away 
from 1 as ¢ increases), while x = —1 is a stable solution (points 
near x = — 1 move towards — | as ¢ increases). 

For a further example involving one variable, consider the 
equation 


This occurs in physics, where x denotes the amount of an element X 
subject to radioactive decay. The equation states that, at any given 
time t, the rate of decay dx/dt is proportional to the amount x(t) 
which remains. Here A is a positive real number, called the decay 
constant. If the amount of X at time t=0 is xo, the equation has 
the unique solution 


x = xge~t. 
This solution can be represented by a single diagram 


0 x. 
mersmsmumanniifes ster uammenvermemusmemg trout 


showing that x decreases to zero as ¢ increases. Notice that since 
by definition x > 0, only the right half of the line is of interest. 
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It is this kind of simplified diagram which is used to represent 
solutions of autonomous systems in two variables. We give an 
example which also arises when considering radioactivity. 

Consider three elements X, Y, Z such that X decays to Y (with 
decay constant A>0) and Y decays to Z (with decay constant 
p> 0). If x(t), y(t) denote the amounts of X, Y which remain at 
time ¢ then 


As before, we represent x, y by a point in R?; notice however 
that since by definition x > 0 and y > 0, only the top right quadrant 
is of interest. 

Now assume A < p. What can be said about the relative values 


of x, y? 
The equations imply that if x 40 
dx —Ax Ax 


Therefore the movement of the point (x, y) can be represented by 
curves with the following properties: 
(i) along the line y = 0 each curve has slope — 1; 
(ii) between the lines y = 0 and y = (A/y)x the slope is negative; 
(iii) between the lines y = (A/u)x and x = 0 the slope is positive. 
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Moreover there is a solution along the line y= kx if and only if 
k = —1+(p/A)k; that is for 


There is also a solution y = yoe~“* along the line x = 0. These facts 
are summarized in the diagram. Notice in particular that if the initial 
position Xo, yo has yo < (A/)xo then y will at first increase and then 
decrease. 


4.2 Phase space 


A large class of examples of autonomous systems in two variables 
arise from second-order differential equations. Let t+>x(t) be a 
function representing the movement of a point along a line, and 
subject to a differential equation of the form 


dt? ae 


where the function / is independent of t. Let y= dx/dt. Then x(t) 
represents the position at time ¢, and y(t) the velocity at time ¢. 
The above equation is equivalent to the autonomous system 


It is now possible to plot solution curves in the (x, y)-plane, and 
to show that, if the function h satisfies suitable conditions, there 
is a unique solution through each point x0, yo. The (x, y)-plane, in 
which the distance x is plotted against the velocity y is often called 
phase space. The equation dx/dt = y implies that as ¢ increases, the 
point (x(t), y(t)) will move ‘clockwise’ along solution curves: from 
left to right in the top half y > 0 of the plane, and from right to left 
in the bottom half y <0. 


4.3 Examples 
Consider the equation 


d*x 
qe + Px =: p>0d 
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which arises from harmonic motion. It can be written as an auto- 
nomous system 


The equations imply that the origin x = y =O is a rest point. In 
order to plot solution curves note that, if y 40, 


and therefore px + y(dy/dx) =0. Along y=0 the slope is infinite. 
By 3.6 there are solution curves px?+y2=c. As t increases, the 
point (x(t), y(t)) moves in a clockwise direction around an ellipse 
px*+y?%=c with centre the origin. The value of the constant c 
may be determined from the position xo, yo when ¢ = 0, since then 
c= px2 + yr. 

An important special case arises when p = 1. Then each ellipse is 
a circle x2 + y2=c. 


_. ax 
“laa” 


harmonic motion 
A second-order equation which arises more frequently in practice 
iS 


d*x dx 
aay + px = 0, 


where p >0 and q#0. It is often called the equation for damped 
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harmonic motion if g>0, or negatively damped harmonic motion if 
q<0, provided in both cases that q is small. It can be written as an 
autonomous system 


AO es 
eae! Ne aa P qy. 


The equations imply that the origin x=y=0 is a rest point. As ¢ 
increases, the point (x(t), y(t)) moves in a clockwise direction and 
since 


d dx d 
Gy (Pe + Y) = 2px a + dy rae ee 


it moves inwards or outwards according as g>0 or q<0. If q?<4p 
the resulting curves are spirals, as is proved in 4.5. 


px+qy=0 


q>090 q<9 
damped negatively damped 
harmonic motion q?<4p harmonic motion 


In all these examples we shall obtain in 4.5 explicit formulae for 
x(t), y(t) in terms of t. These imply that there is a solution curve 
through every point xo, yo; however, only four such curves have 
been drawn on each diagram. 

Often it is not possible to obtain an explicit formula for solutions 
of a differential equation, and in such cases the kind of curve 
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sketching used in the above examples must suffice. For example 
consider the second-order equation 


: 
ga (1-28-(F) | +x=0 


dt? dt dt 
where b > 0. It can be written as an autonomous system 
cong 
dt 
dy 


ER ‘ale eel cla 
a x+b(1 — x2 — y?)y. 
The equations imply that the origin x = y=0 is a rest point. A 
comparison with the equation for harmonic motion shows that 
moreover the unit circle is a solution curve (since if b(1 — x2 — y2) =0 
the system becomes dx/dt=y, dy/dt= —x for which x2+ y2=c 
is always a solution curve). The term q = — b(1 — x? — y?) is negative 
for x, y inside the unit circle, and positive for x, y outside the unit 
circle. Inside the circle x2 + y2 = 1 the behaviour of solution curves 
is similar to that for negatively damped harmonic motion. Outside 
the circle x? + y2=1 the behaviour is similar to that for damped 


Lae 
a ale 
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harmonic motion. The solution x = y=0 is unstable, since points 
near (0, 0) move away. The solution curve x2+y2=1 is stable 
since nearby points spiral in towards it. 


4.4 Linear systems 


An autonomous equation of the type considered in 3.3 is said to 
be linear if it has the form 

o.., 

ae 
for some real number r. We should warn that in 3.7 the word linear 
is used in a slightly different sense, although both have in common 
the fact that no higher powers of x occur. Note that if r, s are real 
numbers an equation of the form dx/dt = rx +s can be reduced to 
the previous form by a change of coordinates which replaces 
x+(s/r) by x. The results of 3.7 show that every solution of 
dx/dt = rx has the form x = xoe"'. The origin x =0 is a rest point. 
If r>0 the movement is away from the origin; if r <0 the move- 
ment is towards the origin. Essentially only three distinct situations 
occur, which are summarized in the diagram. 


r>0 r<0 
unstable rest point stable rest point 


r=0 
line of rest points 


An autonomous system in two variables is said to be linear if it 
has the form 


dx 

ae : 
2 en 

dt px — qy 


where p, q, r, s are real numbers. Note that we no longer assume 
y = (dx/dt) as in 4.2 and 4.3. This would correspond to choosing 
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r=0, s=1, which is in fact the most important special case of (*) 
because it is equivalent to the second-order equation 

d*x dx 

aie 4+ @ en + px =0. 

A solution to (*) will be a pair of functions t+>(x(t), y(t)). In 
accordance with the examples of 4.2 and 4.3, the point (x(t), y(t) 
will be imagined as moving along a solution curve in the (x, y)-plane 
as t increases. Our aim is to list the various patterns of solution 
curves which can arise. For one variable we have just seen there 
are essentially three distinct situations; in the two-variable case 
we list fourteen distinct situations in 4.5. The basic method is 
first to seek straightforward solutions, and then to build up the 
remaining solutions by taking linear combinations of these. The 
following lemmas, therefore, are devoted to results on linear com- 
binations of solutions. 


Lemma. Suppose that xi(t), yi(t) and x(t), ya(t) are two solutions 
of (*), and that ci, c2 are real numbers. Then 


c1xi(t) + coxe(t), — cayi(t) + caye(t) 


is a solution of (*). 


Proof. Let x = c1x1 + cox2 and y=ciyi+ ceye. Then 


ae 1 de dx 


= €1(rx1 + Sy1) + C2(rx2 + Sy2) 
= rx + Sy. 
Similarly dy/dt = — px —qy. Therefore x, y is a solution of (*). 


In the next lemma we consider for the first time complex solutions 
of a differential equation. It is often convenient to consider a pair 
(uo, vo) of real numbers as a single complex number Xo = uo + ivo. 
Similarly a pair of functions t+->(u(t), v(t)) can be regarded as a 
single complex function t+>x(t), t¢Y. Here @ is, as before, a 
subset of R, and x(t) is a complex number. We write x =u + iv 
and dx/dt = du/dt + i(dv/dt), and call the functions re x =u and 
im x =v the ‘real’ and ‘imaginary’ parts of x. A complex function 
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is therefore defined by a description of its real and imaginary parts. 
For example, if x = u + iv is a complex number we define 


et = eutiv = et(cos v +i sin v) = e” cos v + ie” sin v. 


There are historical reasons for the names ‘real’, ‘imaginary’ and 
‘complex’, but essentially all that is involved is a convenient algebraic 
mechanism for handling pairs of numbers and pairs of functions, 
by using the ordinary properties of addition and multiplication 
together with the requirement that i2 = — 1. 

A pair of complex functions x(t), y(t) is said to be a complex 
solution of (*) if, with the above definition of derivative, 


rx ty 
dy 
ea px — qy. 


Since p, g, r, s are real numbers, it is almost immediate that this 
statement is equivalent to the following: both re x(t), re y(t) and 
im x(t), im y(t) are solutions of (*). The previous lemma now 
implies 


Lemma. Suppose that x(t), y(t) is a complex solution of (*), and 
that a, b are real numbers. Then 
are x(t)+ dim x(t), are y(t) + bim y(t) 


is a solution of (*). 


4.5 Classification of linear systems 


As in 4.4 we use the notation (*) for the linear autonomous system 


dx 

— == f +5 
a 
dy 

—$ = —— DX — 
dt Pp qy 


where p, g, r, s are real numbers. The equations imply that the 
origin x = y=0 is a rest point. Other solutions x(t), y(t) may be 
regarded as curves in the (x, y)-plane. The aim is to classify all 


60 ELEMENTARY DIFFERENTIAL EQUATIONS 


patterns which occur, to show how these depend on the numbers 
P; 7, ’, S, and finally to obtain explicit formulae for the solutions. 


The first step is to ask: is there a point x1, yi for which there 
is a solution along the line joining x1, y1 to the origin? Suppose 
for a moment that such a point exists. Then the corresponding 
solution must have the form 

x = xif(t) 
y=yif(t) 
where t+->/(t)is a function. Then dx/dt = xi f’(t) and dy/dt = yif'(t) 
so that substitution in (*) gives 
xif'(t) = (rx1 + syi) S(O), 
yif (t) = (— x1 — 91) f(0). 
It follows that f(t) is some constant multiple of e™ where 
, mx1 = 1x1 + syi, 
myi = — pX1— qj1. 
Eliminating x1, yi we conclude that m must be a root of the equation 
m? + (q—r)m+ (ps —qr) =0. 
We call this the characteristic equation of the linear system (*). 

Conversely, suppose 71 is a root of the characteristic equation. 

Except in the case r—m,=s=0, any x1, y1 such that 
(r—m)x1+ sy1 =0 

automatically also satisfies the further equation px1 + (g + mi)y1 = 0. 
Therefore there is a solution 

x = x1e™1! 

y = yiem! 
of (*). If m1, meg are two distinct roots of the characteristic equation, 
then there is another solution 

X = xgeMel 

y = yeemat 
of (*). Applying the first lemma of 4.4 we have 


Theorem I. Suppose that the characteristic equation has two 
distinct roots m1, m2. Let (x1, y1) and (xe, y2) be non-zero points such 
that (r — m1)x1 + sy1 = px1 + (¢+mi)yi1 =0 


(r — m2)x2 + Sy2 = px2 + (q + ma)y2 = 0. 
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There are solutions 
x = c1x1e™1! + Coxgeme! 
y=ciyre™! + coygeme! 
of (*), where ci, C2 are arbitrary real numbers. 
l(i) m>me>0 | I(ii) O> m1 > me 
unstable node stable node 


XN 


AVA 


wey 


C>=0 


C.= 


I(iii) mi > m2=0 I(iv) 0 =m, > me 
unstable line of rest points stable line of rest points 
X22 @ 1 9 
i — oS 
Sadist, ‘ail aclaelesieh a 
C5> 0 Se te ee c,;=0 D9 eee 
pe EW NERA Gee SUEUR OOS a a X52 
wee 6 ene Gree aemraricromangimneinenaene | BPS — 
paesetarnncamneypye scene) _setgosne re aA IPRe an aaa CEE MEE Gracstnemmonen @igreteecmonen eateete eerie uniilfascenerne tess ED 
c= 0 y 3 ica 0 


I(v) mm >0O> me 
saddle point 


*2.Y2 
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Remark. Since mm, 4me2 the points (x1, y1) and (xe, ye) lie on 
distinct lines through the origin. Therefore given xo, yo it is possible 
to find a unique solution of the above form such that x = xo and 
y=yo when t=O. In fact it can be proved that if m1 4 me every 
solution to (*) is of the above form. The curve patterns in the 
(x, y)-plane corresponding to the above solutions are of five basic 
types I(i)-I(v). Note that if m1 > me, and if c1, ce, m1, m2, X1, Xe are 
all non-zero, 
dy — ciyymye™1t + coyemge™eat 


dx  cyxymjye™1! + coxemoe™at 


_ c1yimye(™1—Ma)t + coyome 


CiX mye(M™1—Mat + C2xX2M2 


Thus dy/dx approaches 


C2yv2MNe2 ey y2 
Cox2Me2 xe 


the slope of the line through (x2, y2), as t becomes large and negative. 
Similarly dy/dx approaches 


Ciyimy Hn yl 
CiX1My xy 


the slope of the line through (x1, y1), as t becomes large and positive. 
We leave the reader to make appropriate adjustments to these 
statements when any of c1, c2, m1, me, X1, X2 are zero, and to check 
that the diagrams I(i)-I(v) can be drawn using the above informa- 
tion. 

In order to complete the classification of solutions of the linear 
system (*) it is necessary to consider the possibility that the charac- 
teristic equation does not have two distinct real roots. In all there 
are three cases. 


Case I: two distinct real roots m, > mp. 
Case II: two distinct complex roots m, m. 
Case III: one repeated real root m. 


Case I occurs if (¢+r)2>4ps, and has been considered above. 
We now turn to case II, which occurs if (¢ +r)? < 4ps. In this case 
therefore s4#0. Suppose that the characteristic equation has a 
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complex root m=f+ig, where g#0. Then the other root must 
be m=f—ig. Let k=(m-—r)/s. Then the complex number k 
can be written in terms of real and imaginary parts, say k =u + iv. 
(In fact u=(f—r)/s and v= g/s.) 
By the same argument as was used in the case when m is real, 
there is a complex solution of (*) 
x=em 
y=kem 
which can be separated into real and imaginary parts 
x = eft cos gt+iest sin gt 
y = eft (ucos gt — v sin gt) + ieft (v cos gt + u sin gt). 
Applying the second lemma of 4.4 we have 
Theorem II. Suppose that the characteristic equation has two 


distinct complex roots m, m. Letm=f + ig andk =(m—r)/s=u+t iv. 
There are solutions 


x = aeft cos gt + best sin gt 
y = aett (ucos gt — v sin gt) + bef (v cos gt + usin gt) 


of (*), where a, b are arbitrary real numbers. 


Ili) f>0 Ili) f< 90 

unstable focus stable focus 

‘ ee 
Ili) f =0 


centre or vortex point 
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Remark. Since m #m we have v #0. Therefore given xo, yo it is 
possible to find a unique solution of the above form such that 
x = Xo and y= yo when t=0. In fact it can be proved that if m4 m 
every solution to (*) is of the above form. The curve patterns in 
the (x, y)-plane corresponding to the above solutions are of three 
basic types II(i)-II(iii) illustrated in the diagram. The flow may 
be in a clockwise or anticlockwise direction; the diagram is drawn 
with flow clockwise. 

It remains to consider case III, in which the characteristic equation 
has a repeated root m. Then (q + r)? = 4ps and m = }(r —q), so that 
also (r —m)? =(q+m)?2=ps. | 

One possibility is that p= s=0 and m=r= —gq. In this case (*) 
becomes 


rs i 
dt 
dy 
~- = WY, 
PEE tad 
Illi) m>0 II(ii) m<0 
unstable node stable node 


ee ae 
AS AS 


UI(iii) m=0 
plane of rest points 
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Then (xoe™!, yoe™*) is a solution of (*) for any real numbers xo, yo. 
The corresponding curve patterns III(i), (ii) are similar to I(i), (ii) 
but every curve is straight. Again, there is a unique curve through 
each point (xo, yo). In case III (iii) the only solutions are constants 
x = x0, ) = Yo. 

The other possibility is that p, s are not both zero. In this case 
the repeated root m yields solutions 


x= xem 
view en 


where the point (x1, y1) is determined up to a scalar multiple by the 
equation | 
(r — m)x1 + syi =0, 


or the essentially equivalent equation 
— pxi—(q+m)yi = 0. 


Now choose some point (x2, y2) which is not a multiple of (x1, y1), 
and define 


a, = (r — m)x2 4+ Sy2, 
by = — pxe —(q + m)y2. 


It follows from the relations (r — m)? = (q+ m)? = ps and 2m=r—gq 
that ai, b; satisfy (r—m)ai+sb,=0 and —pa,—(q¢+m)bi=0. 
Therefore (a1, 51) is some multiple of (x1, y1). Thus there is a real 
number f such that a; = fx: and bi = fy1. Consider the functions 


x = (frit + xe)em, 
y = (fit + yaye™. 


It is easy to check that the above relations imply that x, y is another 
solution of (*). By the first lemma of 4.4 any linear combination of 
the two solutions obtained is also a solution. Therefore if ci, cz are 
arbitrary real numbers the functions 


x = (c1 + c2ft)x1e™ + coxeemt 
y=(c1+ ceft)yie™ + coy2ze™ 


give a solution of (*). Again there is a unique curve through each 
point of the plane. The curve patterns III(iv), (v) are similar to 
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IlI(iv) m> 0 Ill(v) m<0O 
unstable node stable node 
XV an 


c2=0 


Iil(vi) m=0 
boundary layer 


- x) | 


C7=0 c,=0 


I(i), (ii) but now the two straight line solutions have coincided, 
and dy/dx approaches the slope of the line through (x1, y1) both 
as t becomes large positive, and as t becomes large negative. For 
III(vi) see example 2.12: points on opposite sides of the line of rest 
points move in opposite directions. 


Remarks. Not all the fourteen types of linear system I (i)-(v), 
II(i)-(iii), III@)-(vi) are of the same significance from a physical 
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point of view. In most physical situations we expect structural 
stability: that is, a small change in the coefficients p, q, r, S leads to 
a small change in the solution curves. 

Clearly a small change in the equation giving I(iii) will yield 
either situation I(i) or I(v). Therefore I(iii) is not structurally stable; 
neither are I(iv), II(ijii) or III(i)-(vi). However, the fact that a situ- 
ation is unlikely to occur in physics does not prevent its being of 
interest to the mathematician (and even to the physicist — as is 
shown by the example of harmonic motion, of type II(iii), and by 
the example of zero acceleration, of type III(vi)). 


4.6 Second-order linear equations 


The detailed results of 4.5 can be applied to the second-order 
equation 
d*x dx 


ag hE Bee Ys te 


to find solutions t+ x(t). Here p, g are real numbers and the cor- 
responding autonomous system is 


a 
dt ys» 


dy _ 
ol px — qy. 


In the analysis of 4.5 we put r=0, s=1. Then the characteristic 
equation of the system becomes 


m2 +qm+p=0. 


This equation is also called the auxiliary equation of (**). If m 
is a root of the auxiliary equation we take x1=1, yi1=m-—r as 


the non-zero point satisfying (r —m)x1 + sy: =0. Then Theorem I 
of 4.5 implies immediately: 


Theorem I. If g2 > 4p the auxiliary equation of (**) has two distinct 


real roots m= —43q+/(q2—4p) and mz= —43q— V(q2—4p). 
There are solutions 


x = cye™1! + coemat 


of (**), where ci, cz are arbitrary real numbers. 
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The corresponding diagram of solution curves in phase space 
will be of 
type I(i), unstable node, if q?> 4p >0,q <0, 
type I(ii), stable node, if g2 > 4p > 0, g > 0, 
type I(iii), unstable line of rest points, if p = 0, q < 0, 
type I(iv), stable line of rest points, if p =0, g > 0, 
type I(v), saddle point, if p < 0. 
In the same way Theorem II of 4.5 implies 


Theorem II. Jf g2<4p the auxiliary equation of (**) has two 
distinct complex roots f + ig, where f= —4q and g=4+/(4p — q?). 
There are solutions 

x = aelt cos gt + bef sin gt 
of (**), where a, b are arbitrary real numbers. 


The corresponding diagram of solution curves in phase space will 
be of 


type II(i), unstable focus, if g2 < 4p, g < 0, 
type II(ii), stable focus, if g2 < 4p, q > 90, 
type II(iii), centre, if p > 0, q=0. 

The final case is dealt with by 


Theorem III. Jf g2 = 4p the auxiliary equation of (**) has a repeated 
root m= — 4q. There are solutions 
x =(c1 + cat)em 


of (**), where c1, C2 are arbitrary real numbers. 


Proof. We put r=0, s=1 in the analysis of 4.5. Since s40 the 
cases III(i), Ili), I1Gii) cannot arise. We choose x1 = 1, y1 = — 49 
and x2=0, ye=1. It then follows that a1 = 1, b1 = — $q and hence 
that f=1. Therefore the final solution obtained in 4.5 becomes 
x = (c1 + Cet)e™. 

The corresponding diagram of solution curves in phase space will 
be of 


type III(iv), unstable node, if g? = 4p, g <0 
type III(v), stable node, if g? = 4p, q >0 
type III(vi), boundary layer, if p=g = 0. 
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4.7, Examples 


(i) Consider the equation 


d2x dx 


The auxiliary equation m? — 3m +2=0 has roots m=1 and m=2. 
By Theorem I there are solutions 


x = cyet + coe?é 


where C1, C2 are real numbers. 
(ii) Consider the equation 


d*x dx 
Cae +x=0. 


The auxiliary equation m?+m+1=0 has roots m= — } + i(/3/2). 
By Theorem II there are solutions 


x = ae-*t cos (: ¥) + be-* sin (: ¥) 


where a, b are real numbers. 
(iii) Consider the equation 


Oe yee 
dt? dt 
The auxiliary equation m2 —2m+1=0 has a repeated root m= 1. 
By Theorem III there are solutions 
x =(¢c1 + caft)et 


where ¢1, C2 are real numbers. 


Exercises 
Rhee 

dx 

aut Rl 

dt ‘i 


be the equation controlling the radioactive decay of an element 
X, as in 4.1. The half life T of X is defined as the time taken for 
half a given amount of X to decay. Prove that 


T=, log he 
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Consider the autonomous system 


nonneane 

dt 

dy 

re ie 2 
: Ar ey 


Show that there are solution curves lying along the lines x = 0 
and x + y=0. Deduce that the solution curve through (x0, yo) 
1S 

x = xoe! 

y = (x0 + yo)e*! — xoe?. 
Use the substitution t = e” to solve the differential equation 


7a “x i 
What is the pile substitution in the region t <0? 


Write down differential equations with the following solutions: 
(i) x =e (acos (t1/3) t+ b sin (t1/3)), 

(ii) x = ae’ + be-2t, 

(iii) x = (a+ btje-*. 

Find all solutions of each of the following differential equa- 

tions. 


dx a2x dx 


i diab is 0, (il) a a ae Os 
dx sil ae 
Gi) —4 Gp toe 0, @v) 5 +3 =0. 


Find a solution of each of the following differential equations 
with the property: x =0 and dx/dt =1 when t =0. 


d*x dx 
De 3 + 2x =0, 
(i) OF + x= 0, 
ci) +24 x40, 
(iv eg 0. 


dt? dt 


CHAPTER FIVE 


Linear equations 


The word ‘linear’ has already been used in three senses: to describe 
first-order equations of the form (dx/dt) + xg(t)=h(t) in 3.7, to 
describe the autonomous system 


dx 
ona rx + Sy 
dy 
at ae ey 
in 4.5, and to describe the second-order equation 
d2x dx 


in 4.6. In this chapter we seek to justify our use of the word by 
placing these three examples in a single framework. The essential 
feature is that given two solutions of the equation, any linear com- 
bination of them will be a solution either of the given equation, 
or of a closely related equation. This fact leads directly to the 
results often summarized as rules about ‘arbitrary constants’, 
‘complementary functions’, and ‘particular integrals’. 


5.1 Linear operators 


Let FZ be a subset of R which is kept fixed in all that follows, but 
which in practice would often be the whole of R. We consider 
functions f+» x(t) with domain Z. The idea behind linearity is that 
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any linear combination of functions with domain @ is itself a 
function with domain &. If ci, cz are real numbers and x1, xe 
are functions with domain J then c1x1 + ceXxe is the function defined 
for all te D by 


(c1X1 + C2X2)(t) = c1X1(t) + C2X2(t). 


Let V be the collection of all functions tr x(t) with domain 9 
which are infinitely differentiable; that is, the r-th derivative d"x/dt" 
exists for every integer r > 0. If x1, x2 are functions in V then so is 
C1X1 + C2X2 because 

ll (c1x1 + Cox2) =C lied a il 

dt? aan ee 
It is often convenient psychologically to think of x as an element, 
or point, of V rather than as a function ¢+~.x(t). This makes it 
possible to consider functions which associate, to each element 
x EV, an element 7(x) € V. We use the notation T:V->V for such 
a function, but give it a special name to avoid confusion with the 
functions f»x(t) which are elements of V. 


Definition. A function T:W->V is a Jinear operator if, for all 
Ci, Ce R and Als ae & V, 


T(c1X1 + C2X2) = c1T(x1) + coT (x2). 


Examples. (1) Let c be a real number. The function x!cx is a 
linear operator c: VV. 

(ii) Differentiation xt-> (dx/dt) defines a linear operator D: VV. 
Similarly the r-th derivative xt+-~(d'x/dt") defines a linear operator 
Dr:V—->N. 

Gu) If geV the x++(dx/dt)+xg is a_ linear. operator 
(D+ g):V—-V. 

(iv) If p, gq are real numbers then 

PU ss +q oe + px 
dt? dt 
is a linear operator (D2+qD+p):V-—V. 
(v) If ai1,...,an€V then 
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is a linear operator (D"+a,D"-1+ ++++an,):V—V. Here 
ai,..., Qy are elements of V (i.e. functions with domain 9). If 
@i,..., Qm are all real numbers (i.e. constant functions) the linear 
operator is said to be of n-th order with constant coefficients. 

The above examples suggest that certain differential equations 
can be written in the form 7T(x)=0, where 7:V-—>V is a linear 
operator. This leads to the following terminology. If T:WV—V is a 
linear operator, the collection Vo of all x € V such that 7(x) =0 is 
called the solution set, solution space or kernel of T. An element 
x € Vo is called a solution of the equation T(x) = 0. 

The definition of linear operator implies immediately: 


Lemma. Let T:V —V be a linear operator and Vo = {x EV: T(x) = 0}. 
If x1, X2 € Vo then cix1 + cex2 is an element of Wo for any real 
numbers Ci, C2. 


Translated in terms of solutions this lemma becomes: if x1, x2 
are two solutions of the equation 7(x)=0 then cix1+cexe is a 
solution for any real numbers cj, Co. | 

This may be compared with the first lemma of 4.4 which results 
from the fact that the function 


d. d 
(x, yro(S me Me OY, 7 nh pe +a) 


is a linear operator V2--V2. Here V? denotes the set of pairs 
(x, y) with x € V, y € V, and linear combination of pairs defined by 


C1(X1, V1) + Co(X2, yo) = (c1X1 + CoX2, C1y1 + C2y2). 


This explains the use of the word Jinear in 4.5. 
More generally the words /Jinear equation are often applied to an 
equation of the form 


T(x) =h 
where 7:V->V is a linear operator, and h is a fixed element of 
V (i.e. a function t+-+A(t) with domain Y), for example the equation 
oe 
dt 


considered in 3.7 is of this kind, with T the linear operator D + g of 
example (111) above. 


+ xg(t) = A(t) 
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Lemma. Let 7:V->V be a linear operator, and y € V a particular 
solution of the equation T(x) = h. Then every solution of T(x) =h has 
the form yo+y, where yo is a solution of the equation T(x) = 0. 


Proof. The definition of linear operator implies that 
T(vo + y) = Tivo) + TO) =O +h=h 


and therefore yo + y is a solution of 7(x)=h. Similarly if z is an 
arbitrary solution of 7(x)=h then 


T(z — y) = T(z) -—T(y) =h-h=0 


and therefore z — y = yo for some solution yo of T(x) =0; in other 
words, every solution z of 7(x) =h has the form yo + y. 


Corollary. Let x1, ... , Xn € V be functions with the property: every 
solution yo of T(x) =0 has the form yo=c1X1+ *** +CnXn where 
C1,..., Cn are real numbers. If y&V is a particular solution of 
T(x) =h then every solution of T(x) =h has the form 


cixit *** +CnXnty. 


There are traditional names for some of the terms which occur 
in the corollary. The function y is called a particular integral of 
T(x) =h. The expression c1xX1 + *** +CnxXn is called the comple- 
mentary function, and the real numbers ci, .. . , Cn are called arbitrary 
constants. 


Example. Consider the equation 


| d2x dx 
5 Saeed = 
(l+t ) + 27 pie 


It is of the form 7T(x)=h where h is the function t+ 2t and T is 
the linear operator (1 + t?)D2 + 2tD. The equation can be written 


d f dx 

ue —— | =2f, 

4 (a+ F) 
Therefore 


dx 
I 2) =] +4 72 
(+e) +%t*°+C) 
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where ¢, is a real number. It follows that 


dx Cl 


coy te Dp a 
dt 1+? 
x=t+c,tan-1t+ ce 


where C1, C2 are real numbers. Here ¢ is the particular integral, and 
c, tan-! t+ cz is the complementary function. The arbitrary con- 
stants ci, co can be determined as soon as additional information is 
specified about the desired solution. For instance, suppose we seek 
a solution to the above equation for which x =0 and dx/dt=2 
when ¢ = 0. Substitution of these values gives c1 = 1, co = 0. There- 
fore the solution with the required values when t = 0 is 


x=t+tan-¢. 


5.2 Linear equations with constant coefficients 
Consider the equation 


d"x d”-lx 


a Oona Bn + anx = h(t) 


where / is a fixed element of V (i.e. a function t+ A(t) with domain 
QY) and ai, ..., An are real numbers. It can be written in the form 


T(x) =h 


where T:V->V is the linear operator D® + a,D"-1+ +++ +ay of 
example (v) in 5.1. By the second lemma of 5.1 it is possible to find 
all solutions of the equation 7(x) =/h in two steps: 


(i) find all solutions of 7(x) =0, and 
(ii) find a particular solution of 7(x) =h. 


A method for carrying out step (ii) is given in 5.4; however, in 
practice it is often possible to guess a particular solution. We there- 
fore concentrate first on step (i). 

As in 4.6 we introduce the auxiliary equation 


m® + aym™ 1+ +++ +an,=0. 
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Lemma. /f the real number m satisfies the auxiliary equation 
m™ + aym™-1+4 +++ +a, =0 then 


x =cemt 


is a solution of the equation T(x) = 0 for any real number c. 


Proof, If x =ce™ 
T(x) =(D" + a,D"-1 + +++ +ay)cemé 
ent D"(ce™*) + qj D*-\(cemt) ccc an(ce™) 
aatat m(ce™t) + aym"™-1(cemt) ooo aAn(ce™*) 
= (m™ + aym™"1 + +++ + an)cem =0. 


It may happen that the auxiliary equation has a complex root 


m=f+ig. In this case the same argument will imply that any 
multiple of 


emt — eftttgt — eft(cos gt +i sin gt) 


is a complex solution of the equation T(x) =0. The same argument 
used already in the second lemma of 4.4 then implies: 


Lemma. /f the complex number m=f+ig satisfies the auxiliary 
equation m™ + aym™-1+ +++ +an,=0 then 
x =aelst cos gt + beft sin gt 
is a solution of the equation T(x) =0 for any real numbers a, b. 
By taking linear combinations of the solutions given by the pre- 


ceding lemmas, we can build up more complicated solutions of the 
equation 7(x) = 0. The result is: 


Theorem. Let T:V—V be the linear operator 
Oa Pet oP an 


where ai,..., Qn are real numbers. Suppose that the auxiliary 
equation 


m” + aym™-1 + eee +a,=0 
has n distinct roots: 2s complex roots fi + igi, ..., fs + igsandn—2s 
real roots mi, ..., Mn-2s. Consider the functions 


x; = efit cos git l<j<s 
Xs+j = efj¢ sin git l<j<s 


Xostk = e™xt l<k<n-2Ds. 
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Then for any real numbers c1,..., Cn the equation T(x)=0 has a 
solution 
X=C1X1 + °°* +CaXn. 


This theorem contains as a special case the solution to the equation 
(dx/dt) + xa =0 obtained in 3.7, and also cases I and II of the 
solutions to the equation (d2x/dt?) + q(dx/dt) + px =0 obtained in 
4.6. To obtain a generalization of case III we must allow for the 
possibility that the auxiliary equation has some repeated roots. 

It can be checked explicitly that, if m is a root of the auxiliary 
equation of multiplicity r, the functions t*e™! are solutions of the 
equation 7(x) =0 for k=0, 1,..., r—1. This implies: 


Proposition. Suppose that in the above theorem xi;= *** = Xi+r-1 
because of a root m of multiplicity r. Then in the expression 
X=C1Xi + °** +CnXn the terms cixit + °° + Ci+r-1Xi+r-1 may be 
replaced by (ci + cizit + °° + Cétr—1t™-Y)X%. 


Remark. It is also true, although we have not given the proof, that 
every solution of the equation 7(x)=0 is of the form specified in 
the theorem when the auxiliary equation has distinct roots. If there 
are some repeated roots then again every solution is of the form 
specified, provided that the expression c1x1 + +++ + CnXn is modified 
as in the above proposition. 


Examples. (i) Consider the equation 
nO oe + 11 © — 6x = 20 cos y 
This has the form 7(x)=/A where T= D?—6D2+11D-—6 and 
h is the function t++20 cos tf. The auxiliary equation is 
0 = m3 — 6m? + ilm— 6 =(m — 1)(m — 2)(m — 3). 
Therefore every solution of 7(x) =0 has the form 
x = cye' + coe! + cgze%. 


This is the complementary function. To find a particular solution 
it is reasonable, since h(t) = 20 cos f, to try out linear combinations 
of cos ¢t and sin ¢. In fact 


T (cos t)=sint+6cost—11sint—6cost= — 10sin¥¢, 
T (sin t) = —cost+6sint+11cost—6sint=10cost. 
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Therefore x =2 sin ¢ is a particular solution, and every solution is 
of the form 


x = cye’ + coe! + cge3t + 2 sin f¢. 


In practice the arbitrary constants c1, c2, cz are often determined by 
the initial conditions, i.e. the values of x and its derivatives for some 
fixed value of ¢. For instance the solution to the above equation 
may be required which has 


we 0), Phan ae, Hess a6 8) foe 


when ¢ = 0. Substitution of these values gives 


O=c1+c2+ C3 
0=c1 + 2c2 + 3c3 +2 
—6=c¢1 + 4c2 + 9c3 


and hence ci = 2, cp = — 2, c3 =0. Therefore the required solution 
is 


x = 2(e' — e%! + sin f£). 
(ii) Consider the equation 


PO EON aN dune ge 


Here T= D?—2D2+ D—2 and A is the function t:>4e’. The 
auxiliary equation is O=m%—2m2+m—2=(m-— 2)(m2 +1). 
Therefore the complementary function is 


x= ce" + cocost+c3 sin t. 


To find a particular integral it is reasonable, since h(t) = 4et, to try 
out multiples of e. In fact 


T (e*) = et — 2e' + et — 2et = — 2et 


and therefore x = — 2e! is a particular solution. Thus every solution 
is of the form 


x = ce"! + co cos t+ c3 sin t — 2e-. 


where ci, C2, c3 are real numbers which may be determined from 
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the initial conditions. For instance the solution with values x = 0, 
dx/dt =0, d*x/dt? = 1 when t=0 has 


O=c1+c2—2 
0=2c1 + c3—2 
1 = 4c) — c2—2 


and therefore c; = 1, co= 1, c3 = 0. 
(iii) Consider the equation 
d°x, Gex ax 
Ee a ee 
Git aaa at me 
Here T= D?— D?— D+1 and h=2. The auxiliary equation is 
0=m3—m%—m+1=(m+1)\(n—1)2. There are solutions of 
T(x) =0 of the form 


x = cye"' + coet + cet. 


However, according to the above proposition, the last two terms 
may be replaced to obtain the more general expression 


x = ce"! + (co + cat)et. 


Clearly x =2 is a particular solution. Therefore every solution is 
of the form 


x = cie-* + coe’ + cgtet + 2, 


where C1, C2, C3 are real numbers which may be determined from the 
initial conditions. 


5.3 Stability 
Consider the linear equation 
Toh 


where T= D” +a,D"-!+4 +++ +an 1s as in 5.2. If yEeV is a 
particular solution then by the results of 5.1 and 5.2 every other 
solution has the form 


RCL Ce) 8 SA Cae ae 


where ¢1,...,°Cn are real numbers and x,..., Xn are functions 
determined by the auxiliary equation. 
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Under what circumstances will any two particular solutions 
become arbitrarily close for large values of ¢? Since any two parti- 
cular solutions differ by a function of the form c1x1 + *** + CnXn 
it is equivalent to ask: for what equations do the functions 
x1, ..., Xn become arbitrarily small for large values of t? 


Definition. The equation T(x) =0 is stable if every solution tends 
to zero as ¢ tends to infinity. 


To make this definition precise needs a discussion of limits and 
continuity. For the present purposes it is sufficient to know the 
following fact about exponential functions (compare Moss—Roberts, 
Theorem 4.5.9): if k > 0 then, as ¢ tends to infinity, r*e™* tends to 
zero if m < 0 and to infinity if m > 0. 


Theorem. The equation 
dP alt 
dt™ aun* 
is stable if and only if every root of the auxiliary equation 
m® + aymn-1+ eee + aAn=0 


has negative real part. 


+ coe tanx =0 


Proof. By the theorem of 5.2 the roots of the auxiliary equation 
determine functions x1,..., Xn. The equation is stable if each x; 
tends to zero as ¢ tends to infinity. But x; has one of the forms 
tkemst, tkeljt cos gyt, tkefst sin gjt, the power of ¢ being zero unless 
there are repeated roots of the auxiliary equation. In each case the 
function tends to zero, as ¢ tends to infinity, if and only if m; <0 
or f; < 0. 

In the case n = 2 this notion of stable links up with the terminology 
already used in 4.5 and 4.6. The statement that 


a +q + px =0 
is stable corresponds to the statement that the autonomous system 
dx 
di pelt gt 
dy 
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has a stable node or focus at'x = y = 0; that is, solution curves run 
in towards the origin as ¢ tends to infinity. The analysis of 4.6 
shows that the equation (d2x/dt?) + q(dx/dt) + px = 0 is stable if and 
only if p>0, g>0; this condition agrees with that of the above 
theorem because the roots of m2+qm+p=0 have real part — $q 
if g2< 4p, and real part — $q + $\/(q2 — 4p) if q? > 4p. 

The importance of stable equations lies in the fact that often a 
solution of 

T(x) =h 


is required only once a ‘steady state’ has been reached; in other 
words only for large values of t. If T(x) =0 is stable then any two 
particular solutions of T(x)=h will become arbitrarily close for 
large values of t, so that for many purposes it does not matter which 
particular solution is chosen: the behaviour of the solution for large 
values of ¢ is independent of the initial conditions. 

For example consider the equation 


d?x dx 
Since 2 >0 and 3 > 0 the equation 
d*x dx 
We +2 an +3x=0 


is stable. Moreover x = 3 is clearly a particular solution. Therefore 
every solution tends to the value 3 as ¢ tends to infinity; this holds 
whatever the initial conditions. 


5.4 Particular solutions 


We consider second-order equations 


where p, qg are real numbers and / is a function tt>A(t). The aim is 
to find a particular solution. 

In practice the best method is inspired guesswork. However, we 
give an alternative method which is surer but often slower; it can 
be extended to higher order equations, and to the case that p, q 
are non-constant functions. This is the method of variation of 
parameters. 


82 ELEMENTARY DIFFERENTIAL EQUATIONS 


Suppose that the solutions of the equation 


have been found in the form 
X= C1X1 + Coxe 


where ci, cz are arbitrary real numbers. Consider functions of the 
form 
X = 01X1 + Vexe 


where v1, ve are non-constant functions. Then 


dx _ avy dv2 dx X1 dx2 
a= (Ge x1 +S? xe) +(o1 mer =). 


Suppose that 
dv, ee 
Voc: esol 


Then we have 


oa (Gee aoe) (v old gen oa 
at? \dt dt § dt dt dt? dt? 
and hence 
d*x dx dv, dx. dv2 dxe d X1 dx X1 
de ae ae -(44 +e) re (oe eae +px:] 
d?x~» dx2 
+ oS dis +q— BP *+ pre) 
_ 1 dx | dvz dxe 
datidt dt at 
It follows that x will be a solution of 
dx dx 
Pa aun dla 
if v1, va satisfy the equations 
dv; dve 
a tel 


doy dies , doo da _ 
aa a 
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Elimination gives the equations 


at dt at 
and 
dvz {_ dxe 8 a 
we ( 1 ae X2 eh = xXih 


which may be solved to give (dvi/dt), (dv2/dt); and then integrated 
to give functions v1, ve. Of course in a specific example it may not 
be possible to perform the integration explicitly. We give two 
examples in which it is possible. 

(i) Consider the equation 


dies TL, Lael ite 

dt? cos t’ 2 2: 
Every solution of (d#x/dt?) + x =0 has the form c1x1 + cex2 where 
x1 = cos ¢t and x2 =sin ¢. Since 


dxe ee dx} a 9 , en 
Ae Oleg ca (cos t)? + (sin t)? = 1 
the equations for v1, ve become (dvi/dt)= —(sin t/cos t) and 


(dve/dt) = 1. Hence a particular solution is 
x = (log cos t)x1 + tx 
= (log cos t) cost+¢sint 
and every solution has the form 
x=c,cost+c2sint+ (log cos t)cost+sint¢. 


(ii) Consider the equation 


dt? dt 


In this case we compare the method of variation of parameters 
outlined above with the method of inspired guesswork. 


Inspired guesswork 


The equation 
G°x..' ax 
wing) Wasa lls PT ee 2 
i + °° pa 4t+ 
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has a solution x = 2t. Therefore the equation 
2 
—~+4+—-2x=—-4t (1) 


has a solution x = 2¢+ 1. 
The equation 

d*x dx 

priliaies Moacceggean, | UP 

dt? i Me ly 
is likely to have a particular solution ce’; substitution shows that 
this means (1 + 1 — 2)c =1 which is impossible. As a second guess 
try cte’; this time (c+c+ct)+(c+ct)—2ct=1 and therefore 
c= 1/3. Thus x = (1/3)te* is a particular solution of 


d*x dx ie 
The equation 
d*x dx : 
diz oe me 2x = sin ft 


is likely to have a particular solution a cos t+ 5 sin ft; inspection 
shows that this means 


(—a+6-—2a) cos t+(—b—a-— 2b) sin t=sin t. 
Therefore 6 = 3a and a= — 1/10. Thus 


1 i ne 
me i a COS t— 75 sin # 
is a particular solution of 
d*x dx 
— + — 2x =sin t. 
i + , x = sin (3) 


New using the fact that the equation is linear, we can add the 
particular solutions obtained for equations (1), (2), (3) and assert 
that the sum 


1 1 3 e 
x=2t+1+-te 1 cos t i sin 


is a particular solution of the original equation (ii). 


Variation of parameters 
The equation 
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has solutions c1x1 + cox where x; = e~2¢ and xo = et. Now 


dx2 dx} 
X1 —— — X2 —— =e Alet 4 Qe-2tet = 3e-t 
wi ae 
so that the equations for v1, ve become 
dv rae Pe ees 3 2 
sd a ids (ef + sin t — 41) = — e3t — et sin ¢ + 4te2t, 
ihe i oY i 
4a me ‘(ef + sin t— 4t)= 1+ e-* sin t — 4te-t, 
and corresponding integrals are 


301 = — : est + : e*t (— 2 sin t+ cos ¢) + 2te%t — et, 


3vg=t— : e~t (sin t+ cos t) + 4te-# + 4e-t, 


Thus the particular solution obtained for the original equation (ii) is 


1 | |! NN 
ee ale ON a gh Ny 
( 5° +75 ( sin t-+cost)+—t ;) 


1 ae 4 4 
a a ee 
+(5 40 g (in t+ e081) +5 043) 


1 1 1 3 
= aw ot 4 De Oe A 
are titste 19 %8 # - 19 12 ¢ 


Note that the two methods yield particular solutions which differ 
by (1/9)e!; this is a solution of (d2x/dt?) + (dx/dt)—2x=0 in 
accordance with the second lemma of 5.1. 


Exercises 


1 Let T7:V—V be the function which associates to each x € V the 
element 7(x)=y €V described below. In which cases is T a 
linear operator ? 


j) WO =1, Gi) y) =0, 

Git) (0) = (2, (iv) W(t) = °(dx/d0), 
t 

(v) x)= | xodu, (vi) (0) = 1. 


0 
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Prove that every real number c defines a linear operator which 
associates to each x € V the element cx. 

Prove that if S, T are linear operators then so are the functions 
ST and S+T defined by 


(ST)(x) = S(T(x)), 
(S + T)(x) = S(x) + T(x). 


Deduce that 77 =7...T (r factors) is a linear operator, and 
that if a, b are real numbers then 


(T+ a\(T+b)=T?4+ (a+b)T + ab=(T+ b\(T +a). 


Let T be the linear operator which associates to x € V the 
function t+ tx(t). Prove that DT -TD=1. 


Prove that if m is a real number (D — m)(e™+g) = e™ Dg for any 
g € V. Deduce that (D — m)'(e™g) = e™t Dg. Hence show that 
x=(cotcitt+ ++* +¢r-1t? em is a solution of the equation 
(D — m)'x =0 for any real numbers ¢c1,..., Cn. 


Prove that the equation m3 + m?+m-+6 has a root m= — 2. 
Show that the equation 
dx  d*x . dx 
i 6x =0 
dt8 af dt® wie dt ue 


is not stable (even though the auxiliary equation has positive 
coefficients). 


Use the method of variation of parameters to find a particular 
solution to the equation 

pole MEP aA 

dt? (1 + e)2 
Find all solutions of the following equations; also the solution 
with x = 0, dx/dt = 1 when t=0. 


2 
iL abil (i) Ge x= 
dx d*x _ 
lle sui =: et ists = 
Gig 7 ae Me a + 4x =etcos 2t, (iv) WP Ser sin 2t, 
ey C2 4 alee an (vi) Ta + x= Cost 


dt? dt 


CHAPTER SIX 


Guide to further 
reading 


The reader interested primarily in applications of differential equa- 
tions will find his future reading determined by the fields in which 
he works. In chemistry, economics, engineering and physics there 
are many textbooks dealing in greater detail with particular differ- 
ential equations which arise, and with the methods needed for their 
solution. The present chapter is intended rather to give an indication 
of the further development of the mathematical theory, and to 
explain why anyone interested in differential equations should also 
be interested in linear algebra and analysis. 


6.1 Linear algebra 


Many of the results of Chapters 4 and 5 are special cases of theorems 
about vector spaces. For the present discussion it is not necessary 
to know the precise definition of a vector space. The essential 
feature is that, if v and w are two elements of a vector space W, 
then so is the linear combination av+ bw for any real numbers 
a, b. Examples of vector spaces which occur in the previous chapters 
include: 


R set of all real numbers x, 

R2 set of all pairs v = (x, y).of real numbers, 

R” set of all n-ples v =(x1, ..., Xn) of real numbers, 

V_ set of all infinitely differentiable functions ¢++x(t) with 
domain J, 

V2 set of all pairs v = (x, y) with x, y EV. 
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If W, W’ are two vector spaces, a function ¢:W—-W’' is a rule 
which associates an element ¢(w) € W’ to each w € W. The function 
p:W—W' is a linear map if 


$(av + bw) = ag(v) + bd(w) 


for all a, bE R and v, we W. For example the linear operators 
of 5.1 are linear maps T: VV. Similarly the function 


(x, y)t-> (rx + sys) + px. — gy) 


of 4.4 is a linear map ¢:R2->R2. Many theorems about differential 
equations are best understood in terms of results on vector spaces 
and linear maps. For this purpose the following concepts are 
important. 

The kernel of a linear map ¢:W—-W’‘ is the set 


Wo = {w € W:¢4(w) = 0}. 


The definition of linear map implies that ¢(v) = ¢(w) if and only if 
v — wis in Wo. The set Wo is itself a vector space, since if v, w € Wo 
then av + bw € Wo. The lemmas in 5.1 are special cases of these 
properties. 

The dimension of a vector space W, if it exists, is the minimum 
integer n for which the following result holds: there exist elements 
W1,..., Wn € W such that every w € W can be expressed uniquely 
in the form 


W=C1W1+ °°* +CnWn 


where c1,..., Cn are real numbers. If no such integer exists then 
the dimension of W is said to be infinite. For example the dimension 
of R2 is 2, since any w € R2 can be expressed uniquely in the form 
x(1, 0) + (0, 1) where x, y € R; similarly the dimension of R” is n. 
On the other hand the dimension of V is infinite. The main theorem 
of 5.2 states that, if 7: WV —> V is the linear map 


ae LR NR i lag 


where a1, ..., @n © R, then the kernel of T has dimension n. The 
concept of dimension is needed to make precise intuitive ideas 
about ‘arbitrary constants’ and ‘degrees of freedom’. 

The eigenvectors of a linear map ¢:W-—-W are those elements 
w € W such that mw = ¢(w) for some real number m. We give two 
examples to illustrate the fact that many properties of ¢ may be 
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deduced from corresponding properties of the eigenvectors of ¢. 
Firstly, let D:W—>V be the linear map x+~>dx/dt of 5.1. Its eigen- 
vectors are the functions x =ce™! where c, meER. Secondly, let 
¢:R2-—>R2 be the linear map (x, y)t>(rx+ sy, —px—qy) of 4.4 
Its eigenvectors are the elements (x1, y1) € R2, such that 


MX, ='rx1 + Sy1 


ELS PAL or PL 


for some m € R, which played a central role in the.proof of 4.5. 

The basic properties of vector spaces, linear maps, kernels, 
dimension and eigenvectors may be found in many textbooks; for 
example 


K. Hoffmann and R. Kunze, Linear Algebra (Prentice-Hall, New 
Jersey, 1961), 

N. H. Kuiper, Linear Algebra and Geometry (North-Holland, 
Amsterdam, 1962). 


6.2 Analysis 


It has been mentioned in 1.7 that a precise definition of differentiation 
and derivative requires the concept of a continuous function 
f:9->R. Here the domain @ is contained in R, and f is a rule 
which associates a real number f(t) to each t € Y. If D is an interval 
then every differentiable function is continuous, and every con- 
tinuous function is integrable. However, the converse statements 
are false. All properties of differentiable functions, including assump- 
tions 1.10 and 1.12, the chain rule, the rules for differentiating 
sums and products, and the precise definitions of trigonometric 
and exponential functions, depend ultimately on the properties of 
continuous functions developed in a first-year university analysis 
course. Suitable references are 


R. M. F. Moss and G. T. Roberts, A Preliminary Course in Analysis 
(Chapman and Hall, London, 1968), 

T. M. Flett, Mathematical Analysis (McGraw-Hill, London, 
1966). 


The latter book discusses in addition the continuity and differen- 
tiability of functions f:4->R” where the domain F is contained 
in R™, This leads to the definition of partial derivatives, and to such 
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results as the chain rule and inverse function theorem for functions 
of several variables. The study of such functions leads to partial 
differential equations; these have a vast literature and many applica- 
tions, but a discussion of these would lead far away from the scope 
of this book. We mention, however, three examples in which func- 
tions of this type have already occurred in earlier chapters. 

(i) The differential equation 


* f(t,x), (Lne”, 
dt 
considered in 3.5 involves a function /:e7? +R where e” is a region 


in R2, This is the case m= 2,n=1. 
(ii) The linear system 


dx 

— =X S 

Fe ited 
ay 
eee 


considered in 4.4 can be written in the form 
2 = $0), v=(%, WERY, 


where ¢:R2->R2 is the linear map mentioned at the end of 6.1. 
Solutions are functions v:R-—>R? which associate to each teR 
the point v(t) = (x(t), y(t) € R?. 

(iii) The differential equation 


d™x . d"-1x 
POEL EEC 1 POS CSU ek EET a Se 
dt” dtn-1 


can be turned into a system of first-order equations by substituting 
X1= x, Xe=dx/dt,..., Xn =d"-1x/dt*-!. It can then be written 


+ +++ +anx =0 


dx1 « 
eee ae ~~ 2 
dt 
dxXn-1 
dt i 
Xn 
—— = — AXn- °** —aAnX1. 
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If 6:R”—R” is the linear map which associates to 
== (xi, ie » Xn) E R® 
the element 
P(v) = (Xe, bon gi the To aLAg =." '* * — AnX1) 


of R”, the equation T(x) = 0 of 5.2 is equivalent to the equation 
Sade, | v= (tu...) xe ER 


Solutions are functions v:R->R®” which associate to each t ER the 
point v(t) = (x1(t), ... , Xn(t)) € R*. 


6.3 Existence theorems 


The main question left open in Chapter 3 is: for what regions co” 
in R2, and for what functions f:c” R, does the equation 


dx 
dt = f(t, x), (t, x) E oS, 


have a solution through every point (to, xo) € &”, and when is this 
solution unique? Theorems 3.1 and 3.3 answer the question in the 
special cases f(t, x) = g(t) and f(t, x) = A(x). What is required is a 
theorem which does not give such complete information in these 
special cases, but which on the other hand applies to a more general 
class of functions f For simplicity we assume that e” is a rectangle 
{(t, x) € R8:p<t<q and a<x<b} in R?. 


Existence theorem. Let f:c’->R be a continuous function, and 
suppose (to, Xo) € &%. There is a real number h>0 such that the equa- 
tion 


dx 
dt = f(t, x) 


has a solution t'>x(t) with domain G={tER:ty—h<t<tot+h} 
and x(to) = Xo. 


Uniqueness theorem. Let f:07 —R be a continuous function for which 
the partial derivative 0f/0x is also continuous, and suppose (to, Xo) € ©”. 
If the equation 


dx 
=f, %) 
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has two solutions x and x with x(to)=X(to)=x0 and domain 
BG={tEeR:ty—h<t<to+h} then x(t) = X(t) for allt eZ. 

These two theorems imply immediately the truth of assumptions 
1.10 and 1.12. The necessity of the additional hypothesis (con- 
tinuity of the partial derivative 0f/0x) in the uniqueness theorem is 
shown by examples 2.5 and 3.4. 

Both theorems have immediate extensions to the case of a 
differential equation 


dv 
dt = f(t, v) 


where 0 = (X1,...5 Xa) ER and ff, 0) =(1, x1...) Xa ee ae 
region ec? of R”+!, It is this extension, applied to example (ii) of 
6.2, which shows that the linear system of 4.4 has only those solu- 
tions found in 4.5. The same extension, applied to example (iii) 
of 6.2 shows that every solution of the linear equation with constant 
coefficients 7(x) = 0 is of the form claimed in 5.2. 

Other extensions of the existence and uniqueness theorems deal 
with the question whether the solution through (to, xo) depends 
continuously on (to, xo), and with the properties of approximate 
solutions. Here the theoretical discussion is also of great practical 
importance, because it is often possible to compute a numerical 
solution by obtaining closer and closer approximations to it. All 
these matters are discussed in more advanced textbooks on ordinary 
differential equations, such as 


E. A. Coddington and N. Levinson, Theory of ordinary differential 
equations (McGraw-Hill, New York, 1955), 

S. Lefschetz, Differential equations: geometric theory (Inter- 
science, New York; 2nd edition 1962), 

W. Hurewicz, Lectures on ordinary differential equations (M.1.T. 
Press, Cambridge, Mass., 1958). : 


The latter book is also the best introduction to the qualitative 
theory of differential equations, in which one adopts the approach 
of 4.5 and asks (in higher dimensions, or for a non-linear differential 
equation) what patterns of curves can arise, and which equations are 
structurally stable. 


Notes on 
exercises 


Introduction 


1. y can have arbitrarily small values when x <0. 2. A partially 
correct statement is y=sin (a+ x) for —407<a+x < 4m (see 3.4 
for a complete statement). 3. Function uv = 1/x is not defined when 
x =0. 4. First function never makes sense, second function only 
when x = 1. 5. Consider regions x > 0 and x < 0 separately. 6. Solu- 
tion has distinct formulae for x > 0 and x <0. 


Chapter One 


1. Domains of functions given are x >0, x >0, all x, x >0, x >1, 
2nt < x < (2n+ 1)z, empty set. 2. (i) a + 2./x forx >0.(li)a + logx 
for x >0, and b+log (—x) for x <0. (iii) a+sin x for all x. 
(iv) an + log sin x for 2n7m < x < (2n + 1)z, and by, + log (— sin x) for 
(2n — 1)a < x < 2nz. (v) a+ log log x for x > 1 and b + log (— log x) 
for 0 <x < 1. (vi) a+ log (— log sin x) for 2na7 < x < (2n + 1)z. 


Chapter Two 


1. First sketch slope lines along x =0, t=0, x= +1, then fill in 
solution curves. 2. Equations are defined in regions t>0, all ¢, 
t>0, t>0, all ¢, but first two equations are of second degree with 
two solution curves through each point (t, x) with t>0. 3. If t40 
slopes at (t, x) are (1/2t)[x + \/(x?—41)]. Hence no solutions 
through (t, x) if x2 <4f, two solutions if x2 > 4r. If x2 = 4r slope 
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is x/2t = 2/x, that is tangent to parabola x? = 4t. Therefore x* = 41 
is a solution curve. Differentiation of equation gives 

dx\ 2 dx d*x (dx\? d2x 

ent FF AR IAT A iy ER a i 

(5) ' dt dt? (7) [ dt? 

Therefore either 2t(dx/dt) = x (in which case x? = 41) or d?x/dt? =0 
(in which case solution curves are straight lines). Therefore all other 
solution curves are tangents to the parabola. 


Chapter Three 
1. G@) x=a+02+4+ 22? with a= —3 for solution through (1, 0). 
(ii) x=a+t bet with a= — be. (ili) x=a+tan-! ¢t with a= -— fr. 
(iv) x=a+sinh—! ¢ with a= —log(1+ 7/2). (v) x=a+sin ¢ with 
a=-—sin l. (vi) x=an+log sin ¢ for 2n7<t<(2n+1)z, and 
x = bn + log (— sin t) for (Qn — 1)a < t < 2nz with ap = — log sin 1, 
b,, arbitrary. 2. (i) x = 1 or x = 2 or x = (aet — 2)/(aet — 1) wherea 40 
with domain R. (ii) x = — log (a— bt) with domain {t € R:a> bt}; 
choose a =e! to get unique solution through (0, 1). (il) x =1 or 
x=1-(1/t—c) with domain tc. (iv) x=1 or x=cosh (t—c) 
with domain ¢>c; solution through any point on x=1 1s not 
unique. (v) x = 0 or x = (t—c)? with domain t > c; choose c= — 1 
to get unique solution through (0, 1). (vi) x =0 or x =sin~! (ae*) 
where a #0 with domain ¢ < —log | a|; choose a=sin | to get 
unique solution through (0, 1). 3. (i) Either t=0 or x =O or 

1—3t 

t2 


dx 
os | tia 
(2x + x?) e 


which has variables separate. (ii) Since t > 2 either x = — $ or 


Bite ( 1 dx 
12-4 \14+2x/) dt 
(iii) Substitute u = x/t. (iv) Since t > x > 0 equation can be written 


oe 

dt l1—wu? 
where u = x/t. (v) Equation becomes (d/dt)(xe*) = 2t. (vi) Equation 
becomes (d/dt)(t? + 3tx — 4x?) = t?. 4. (i) Write in form 


8 (xe2t) = et, 
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(ii) (d/dt)(x/cos t) =0. (iii) (d/dt)(x cos t) =2 sin 2t. (iv) 


a nee Oe Mee) e-1/t*. 
dt \t8 t3 


(v) Substitute y=x-3. (vi) Substitute v= log ¢. 5. If t= —e¥ 
then 


——— re —er 


dx dt dx _ dx 
Gi dudt’ de. 


6. Existence of derivatives at t=O implies (i) no restrictions, (ii) 
p=r, (ill) p=randq=s. 


Chapter Four 


1. x=Xoe~*t has x = x9 when — At=log $= — log 2. 2. (0, e2¢) 
is a solution along x = 0, and (e#, — et) is a solution along x + y=0. 
Hence (aet, be — aet) is a solution. Choose a= xo, b=x0+ yo. 
3. Substitution ¢ = e” gives 


d2x dx 


and hence x=ae“+be2™=at+ bt? for t>0. In region t<0 
substitute t=—e”. 4. (i) (d2x/dt?) —2(dx/dt)+4x=0. (ii) 
(d2x/dt?) + (dx/dt) — 2x =0. (iii) (d?x/dt?) + 2(dx/dt) + x =0. 5. (i) 
x = ae*t + bet, (ii) x =(a+ bt)e**. (iii) x =e (a cos t +b sin 1). 
(iv) x=a+ be-*. 6. (i) x= —et+e%, (ii) x =sin ¢. (iii) x = te-*. 
(iv) x = et sin ¢. 


Chapter Five 

1. (i) (iii) (vi) non-linear, (ii) (iv) (v) linear. 

2. (ST )(c1x1 + C2x2) = S(T (cix1 + c2x2)) 
= S(c1T (x1) + c2T(x2)) 
= €1ST (x1) + c2ST (x2); 


S + T similar; 
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(T + aT + b)x = (T + a)(Tx + bx) 
=(T+a)Tx + b(T + a)x 
= (T* + aT + bT + ab)x. 
3. DT(x) = D(tx) =x +tD(x)=x+TD(x). 


4, D(e™‘g) = D(e™)g + emt D(g) = e™ mg + e™t D(g). If g is a solution 
(co + cit+ °°* + Cr-it'-!) of D(g)=0 then x =e™g is a solution 
of (D — m)"(x) = 0. 


5. m3 +m2+m+6=(m+ 2)(m? —m-4 3) and roots of m2?—m-+3 
have positive real part. 


6. Choose x1 = et and x2 = et. Equations for v1, ve become 


ik POON AE 052, a 
a ae we On eet 


Integration gives vi =1/(1+ e*). The substitution u=e-¢ yields 


ve = 2 log (1 + e-*) + 1/(1 + e~*) — et. Therefore there is a particular 
solution: 


; si Oy 

x = 2e' log (1 +e) 4+ ~ — wen m9 “ee 

7. (i) x = aet + be-t — t2 —2 where a, b are arbitrary real numbers; 

for solution such that x =0, dx/dt=1 when t=0 choose a= 3/2, 
= 1/2. (ii) x = aet + bet + det with a=1/4, b = — 1/4. (iii) 

x = ae~* cos (t4/3) + be-* sin (t4/3) + (8/73)e! sin 2t + (3/73)et cos 2t 

with a = — 3/73, b =174/3/73. (iv) 


x = ae* sin 1%) + be*t cos (: we) +5 +— cos 2t— hi sin 2t 
13 13 
with a =40/13 1/3, b = — 2/13. (v) 
x = aet + be-3t + i (¢ sin ft—2t cos t) + i. sin t + ice cos t 
10 50 50 


with a= 1/4,b = — 47/100. (vi) x =acos t+b sin t+ 4¢ sin ¢ with 
a=QQ; b= |. 
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line, 7 saddle point, 25, 61 
line of rest points, 57 second degree, 28 
linear combinations, 58, 87 second derivative, 29 
linear equations, first-order, 43 second-order equation, 29, 67 
—, second-order, 67 sin, sin-1, 37 
—., higher-order, 73, 75 sinh, 33 
linear function, 11 slope, 11 
linear map, 88 solution, 13, 73 
linear operator, 72 solution curve, 21 
linear system, 57 solution set, 73 
logarithm, 16 solution space, 73 
square root, 5 
marginal cost, 45 stable equation, 80 
multiplicity, 77 stable rest point, 57 
stable solution, 51, 57 
negatively damped, 55 structurally stable, 67 


node, 61, 64, 66 


tangent, 12 


partial derivative, 41, 89 typical solution, 21 


particular integral, 74 
particular solution, 74 


phase space, 53 unique solution, 17, 44 
plane, 7 uniqueness theorem, 91 
plane of rest points, 64 unstable rest point, 57 


unstable solution, 51, 57 
radioactivity, 51 


real number, 7 variation of parameters, 81, 84 
real part, 58 variables separable, 39 
repeated roots, 64, 68 vector space, 87 


rest point, 50 vortex point, 63 


